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ABSTRACT 


Electromagnetic interference from power lines 1S one of 
the main sources of man-made interference to communications 
in the high frequency (HF) radio Sorae Two types of radio 
interference generated By power lines are gap-type noise 
caused by electric discharges across line hardware and 
corona noise caused by the partial breakdown of the air due 
to the high electric fields around transmission line 
conductors. Using original data, this research has devel- 
oped a parameter based model of gap-type and corona noise 
that allows the fundamental noise mechanisms to be 
mathematically or physically simulated. An expression for 
the power spectral density (PSDB) of gap-type noise and 
corona 1S derived. The energy detection problem is 
formulated, and using analytical results based on the Hall 
model for radio noise, a robust energy detection receiver is 
developed. Tests of this receiver using actual and 


Simulated data are described. 
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I. INTROOUCTION 


A. PURPOSE 

‘Observations of radio interference at high frequency 
(HF) receiver sites have indicated that existing models of 
man-made radio noise are inadequate to describe the observed 
time-. and frequency-domain behavior of the noise. The data 
also shows that quite often a specific noise source pre- 
dominates in a particular location. These observations led 
to the objectives of this dissertation: 


1. to develop specific models for certain man-made 
noise sources in the HF radio band and 


2. to apply this model to the analysis of energy 
detection receivers. 

B. BACKGROUND 

Studies conducted over the past few years have indicated 
that one of the primary sources of man-made radio noise in 
the HF band are alternating current (AC) transmission and 
distribution lines (CRefs. 1,2]. Two of the primary sources 
Of power line noise are gap noise, also known as-~= micro- 
sparking, and corona. Gap noise is caused By a sparking 
process between hardware points on utility poles and corona 
is caused by the partial breakdown of air due to high 
electric fields around high voltage conductors. Both of 


these noise types are non-Gaussian noise processes. 


The majority of the models of man-made and atmospheric 
noise to date Rave assumed the high amplitude impulses 
driving the receiver have either a Poisson arrival rate or a 
variation of Poisson arrivals. This assumption has allowed 
the derivation of first order envelope statistics: the 
amplitude probability distribution (CAPD) Function and the 
density function of the phase of received noise for various 
amplitude distributions of the driving impulses. Observa- 
tions of man-made noise from power lines, however, has 
indicated that the assumption of interpulse independence 
(Poisson arrival times) is not valid. There exists a 
definite time domain correlation of the impulses driving the 
receiver. This is due to the underlying deterministic 
mechanism of the fundamental frequency of the power line 
voltage. Consequently, this research has concentrated on 
the statistics of the impulse arrivals which are manifest in 
the autocorrelation function and the spectrum of the 
observed noise process. Accordingly, the models developed 
are specific rather than general and are more suited to 
source identification by spectrum analysis and robustness 
evaluation of systems rather than to generic specification 
of optimum receiver structures. 

The performance of energy detection receivers’ for 
stochastic signals in non-Gaussian noise has received 
relatively little study. Actual implementation of algo- 


rithms for signal detection has been based on heuristics. 
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The advent of digital signal processing techniques that 
allow for complex post-detection algorithms suggests a 
careful statistical analysis of the noise performance of 
these systems may help develop more efficient and robust 


receivers. 


C. CONTENT 

A brief outline of the rest of the thesis will now be 
presented. Chapter II develops a general model for bandpass 
impulsive phenomena using a filtered impulse model and 
complex envelope theory. Well known empirical models and 
physical models of radio noise are presented and their 
relationship to the generalized model discussed. Particular 
attention 1s paid to the Hall model for atmospheric radio 
noise CRef. 3]. Some useful extensions to it are derived 
that will be used for simulations in Chapter V. 

In Chapter III, field observations from sources of gap 
noise interference are presented and analyzed. Three cases 
of actual interference are used and a probabilistic model 
based ona statistical analysis of the data and the filtered 
impulse model developed in Chapter II is specified. The 
power spectral density (PSD) of gap type noise as predicted 
by the model is derived and compared to the data. 

Using an actual case of corona noise from a 500 kV power 
line, the same type of analysis is carried out for corona 


noise. The corona noise model is based on the same filtered 


1@ 


impulse framework as the gap noise model of Chapter III 
however, the specifications differ substantially. 

Using the Hall model for atmospheric noise discussed in 
Chapter II, a locally optimum receiver for detection of 
unknown Signals in HF atmospheric noise 1S derived. A 
practical modification to this receiver shows it to be an 
adaptive limiter. The performance of this receiver is then 
examined in Simulated man-made noise and in recorded HF 
Signal and noise data. 

Appendix A presents a brief description of the 
instrumentation that was used to collect much of the data 
for this thesis. The dissertation concludes with a brief 


summary of results and some suggestions for future research. 


1 | 


II. IMPULSIVE NOISE AND RECEIVER MODEL 


A. INTRODUCTION 

In this chapter a general set of specifications for an 
impulsive noise interference process will be developed and 
the results of previous work in the field explained in terms 
of the specifications. This generalized scenario will be 
used in later chapters to specify complicated sources of 
man-made radio interference. In attempting to describe 
atmospheric radio noise, two general types of models have 
been developed: empirical models designed to fit first order 
statistical data, and physical models directly related to 
the underlying physical mechanisms. The Hall empirical 
model for atmospheric radio noise has been shown to fit 
atmospheric noise data very well and will be used in this 
dissertation. Some extensions to the Hall model will be 


developed in this chapter. 


B. GENERALIZED SCENARIO 
In order to provide a framework for the discussion of 
the impulsive noise models, a general interference scenario 
for impulsive noise will be described (CRef. 4,5]. A typical 
interference scenario consists of the following elements: 
1. a source of interference, 


Cc. a transmission medium to the receiver and 


Le 


3. the receiver where the interference manifests 
itself. 


For this study the noise process at the | -aceives: is 
modeled as the sum of a high density (Cin time), low amp- 
litude Gaussian component and a low density, high amplitude 
impulsive term. 

These .elements are shown below, where e’(t) in Fig. 1 is 
the impulsive interference, L’(w) is the frequency response 
of the transmission médium, z’(t) i165 white Gaussian noise 
and H’ Cw) is the combined response of the RF and IF filters 


of the receiver. 


2 Cees 





Figure 1. Interference Scenario 


ee 


The source of impulsive noise iS specified by the 

equation 

N¢t) 

e’Ct) = Des (t-t;) (2-19) 

i=1 
where e6, is the real amplitude of the ith pulse and NC(t) 15 
a unit counting process that generates arrival times, the 
t, ‘s. Using this representation, the impulsive interference 
source 1S described by the probability density function 
(PDF) of each @:, given by pi Ce), and the impulse arrival 
times generated by the unit counting process, N(t). 

In order to simplify later analysis the interference 
scenario will be expressed in terms of its complex envelope. 
The interference process at the receiver detector is n’(t) 
and will be replacéd by its complex envelope equivalent, 


n(t), where 


jwot 
n’(t) = Re€Cn(tde J (222) 
and 
nCt) = ne (6) + jms (tO). (2-3) 
jwot 
The term n(toe is known as the analytic signal 


representation of n’C(t), and nctt) 1s the complex envelope 


representation of n’(t). The reference frequency for the 
complex envelope representation 1S We . When n’C(t) 18 a 
bandpass process, Me (t) and met) are the lowpass inphase 
and quadrature modulation components of n(t) respect- 


ively.CRef. 6:p. 75] 
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Specifying a filtered impulse model in complex envelope 
form requires deriving an expression for the driving impulse 
Function given in Eqn. 2-1. Two equivalent approaches may 
be taken and they are both outlined below. In the first 
approach a bandpass impulse may be postulated where the 
spectrum of the impulse is assumed flat from wo - W to Wo + 
W where wo is the reference frequency and W is an arbitrary 
bandwidth. For this case the analytic signal representation 
is easily found in the frequency domain. Using the 
frequency domain definition, the analytic signal associated 
with e’(t) in Eqn. 2-1 is 

se) 
jwt 
e, (t) = ayn | E’ (woe dw. (2-4) 
@ 
Inserting the definition of the bandpass impulse defined 
above and using Eqn. e-1 
NCt) jWeo (t-t, ) 
e,(t) = © Ca/n) e Simcw Cent, ) 17 Cc 6, ). (2-5) 
i1=1 


The complex envelope of e’C(t) is defined as 


-jWot 
e(t) = eg, (te (2-6) 
and 
N¢t) - 39: 
e(t) = £ (2/n) e SintWCt-t, )/CE-ti D C2=7> 
i=1 
where Wo ti = 8.. As W becomes large, e(t) will asympto- 


tically approach CRef. 7] 


iS 


NCt) -j8.1 
e(t) = xX 2e 6:Ce— tee (2-89 
i=1 
The second approach (Ref. 8] is to consider the impulse 
as ideal and constant in the frequency domain. The analytic 
signal representation of e’(t) is then directly determined 


from the time domain definition of an analytic signal. The 


analytic signal is defined as 


e,(t) = e’(t) + je’ (bt) (2-9) 
where @’(t) is the Hilbert transform of e’(t). The analytic 
Signal representation of Eqn. e2-1 using the above defin- 
itions 1S 

N¢t) 
en. (t) = © &CtE-ti) + j/ On (t-ti 1) (2-110) 
i=1 
since 1/nt is the Hilbert transform of the delta function. 
The complex envelope of e’(t) is then 
NCt) = 40% — jWet 
e(t) =X &(t-tide + je /(nCt-t, )) (2=iie 
i=1 
using Eqn. 2-6. 


Although apparently different, when the complex impulse 


trains described by Eqns. 2-8 and 2-11 are convolved with a 


complex filter impulse response, the filter output is 
identical for both representations. If h(t) is the complex 
envelope of the filter impulse response, then the filter 


output with the input given by Eqn. 2-8 is 
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NCt) -j50, 
roe Rce-—t, dD (2-12) 
a 


where complex convolution 1s defined as 


latte eCt) = (172) nica ye Cl—7) da. Caq a. 
ae 
The filter output with the input defined by Eqn. 2-11 15 
NCt) e len = {Wo & 
C172) Rah Ct -caee + je ac ty (2-14) 
1=1 
where the convolution definition of the Hilbert transform 
was used CRef. 6:p. 69]. Using Eqn. 2-6, the same result as 
Mermrved in Eqn. 2-le i185 then found. Due to its simplicity 
and for ease of use in programming simulations, EG. .e=s 
will be used to represent the driving impulse function for 
the rest of this dissertation. 

One further point to mention concerns the probability 
PeeseriDUuUtION OF G6. in Eqn. e2e-S, For the processes we con- 
Sider, t:, will have random arrival times and, as can be seen 
ion eOn. 227, CO \ 1S formed by multiplying ti, times the 
reference frequency, Wo - Consequently, 8, 1s the phase of 
the impulse time with respect to the reference frequency. 
It can be assumed to be uniformly distributed over @ to en 
when the t, ‘Ss have a probabilistic interarrival distribution 


and Wo 1S much greater than the inverse of the interarrival 


times CRef. 9: pp. 279, 10]. 


1 


In an actua ! interference Situation the received 


impulses are filtered by the radiating antenna, spreading 


losses, atmospheric attenuation, recelving antenna and 
cabling losses. These terms will be lumped together as an 
equivalent filter, Ea Goor Examples of this filtering term 


for various interference sources Rave been recorded Dy many 
researchers (CRef. le bead The filter characteristic 1s 
specified in the frequency domain as L’Cw) since that 15 
where 1t 1S most easily observed. 

An example of this type of filtering 1S shown in Fiq. 2 
where the frequency spectrum of an impulsive noise source 
due to agap discharge 1S shown from @ to 200 MHz. The 
Straight lines in the 3-Ax1s view are due to stations in the 
high frequency CHF) band from 2@ to 30 MHz and the FM band 
from 88 to 106 MHz. The continuous envelope seen in the 
upper view 1S the magnitude of the frequency response of 
Ee es The solid line in the upper picture 1s the noise 
Floor of the spectrum analyzer at -10@ dBm. 

The complex envelope of the impulse response of the 
attenuation term, L’(w), 165 given in the frequency domain by 

btw) = Pe be’ <-@ sus) + aera ane (2- ae 
where Le denotes the lowpass part of the quantity ia 
brackets. 

A complex Gaussian component z(t) will also exist in the 


observed noise process and 1s due to the combination of: 
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1. the thermal noise in the receiver, and 


2. the combined sum of many low level atmospheric or 
man-made impulsive sources. 


This term will be modeled by independent zero mean Gaussian 
noise processes for the inphase and quadrature terms with a 
power spectral density (PSD) in each component of No/2. The 
distinction between the Gaussian noise term and impulsive 
noise term is that many receiver responses overlap for the 
high density, low amplitude Gaussian Case, satisfying the 
condition for the Central Limit Theorem. For the impulsive 
noise the receiver response to each impulse 1s discernible 
such that the probability of more than one or two pulses 
overlapping 1S negligible. Common sources of impulsive 
noise are man-made noise due to power lines and atmospheric 
noise due to lightning. 

The time domain input to the receiver 1S now modeled By 
the expression 

ee) 


MiGrD a ONG ween es) 0, 
2 


n(t) = (hCt)/2) = € zee eee ee eCo) aad 


i=1 
oD (2-= Tew 


where 
HCE) = Me Geo + jj neeGe. 


The noise power bandwidth of Hw) will be designated Bere. 
In most situations L’ Cw) will effectively be constant 1n 
comparison to the narrowband receiver filter Hw). Ths 
makes L¢w) a constant that depends only on Wo - 


Additionally, q(t) will be defined by 


2@ 


Ge ecoe— CnCe7e) ) e. 2 Ct) aL, 
and 1s the low level Gaussian noise at the detector. The 
power in the inphase and quadrature components of the 
Gaussian noise will be NoBerr /2. The complex lowpass noise 


process can now be written as 


N(t) = Hele 
n(t) = g(t) + Cae hCt-t, ) (2-18) 
i=1 
where 
a, = e,l’lwo). (2-19) 


The parameter a, is then a scenario dependent parameter that 
is a function of the reference frequency and the impulse 
amplitude. With the above simplifications, Fig. 1 can now 


be expressed in complex envelope form as shown in Fig. 3. 


z(t) 


a(t) H Cw) moe. 


Figure 3. Complex Envelope of Interference Scenario 


The terms from Fig. 1 have Been replaced by their equivalent 
complex envelope representations and the prime is dropped to 


indicate this. 


au 


Some further simplifications) will now be considered 
based on ae characteristics of the H’(w), the bandpass 
fr lbter. If the impulse response H(t) of the filter 15 real, 
H’€w) haS conjugate symmetry about the origin in the sense 
that H’* (w) = H’C-w). Furthermore, if H’(€w) 15 sym-metrical 
in the same way about the reference frequency, Wo 5 then 
Ns (t) will equal @. 

The envelope squared of the process defined in Eqn. 2-18 
1s 


2 2 2 
E (tt) = me (tt) + me (C(t). (2-20) 


If the Gaussian noise term is assumed to equal @ then 


ren NOT) NCCT) 
E <£)-=7 LC 6a A; Me CC=—C, Ihe (tb, cost G, oes 

1=1 j=l 

NOT) NCT) 


is KC a. a, hs (CE—-t, Ds (E-€; 160S(8, -8,). Cee 
1=1 j=l 


If NCCT) defines a low density counting process such that 
there 1s negligible overlap between subsequent pulses and 
htt) is the complex envelope of the impulse response of a 
bandpass symmetrical Filter, then the envelope of n(t) can 
be simply written as 

N¢t) . 

E¢t) = 2 4. Gece (2-ae" 

1=1 

where ACt) 15 equal to he C(t). 


At this point of the analysis two avenues may be 


pursued. If there 1S a time dependency between pulses (non- 


oe 


exponentially distributed interarrival times), then the 
correlation function and spectrum of the process can be 
examined. If the amplitude probability distribution (CAPD) 
of the envelope is the desired result, then the  char- 
acteristic function of the envelope is most useful. The APD 
function 1S commonly used In radio noise research and 15 
equal to one minus the cumulative distribution function 
Cepr). In general these two approaches are mutually 
exclusive since the assumption of time dependency between 
the pulses makes the envelope APD function difficult toa 
calculate. lf the impulse arrival times form a Poisson 
point process then the envelope APD may be calculated: 
however, the spectrum 1S constant. The remainder ofr this 
chapter will Follow the envelope APD approach. Results 
obtained by previous researchers will be developed using the 
general noise model. Chapters III and IV will explain the 
other avenue and look at the spectrum of complicated man- 


made noise processes. 


Gee EMPIRICAL MODELS 

Spi rica i models of atmospheric noise have been 
developed to provide a mathematical expression for the 
first-order statistics of the envelope of the received 
waveform. In particular, the APD has been emphasized. This 
type of model attempts to construct a mathematical express- 


10n that fits observed data without regard for the physics 


ae 


of the interference scenario. Empirical models have the 
advantage that the resulting expressions are much simpler 
than those obtained from the filtered impulse models. One 
disadvantage of this class of model is that since only the 
first order statistics of the noise are considered in 
developing and fitting the model, the RAigher order sta- 
tistics may not match the data well. 

Variations of the Rayleigh distribution have been 
proposed by a number of authors to fit observed atmospheric 
noise data. The Rayleigh PDF is 

=O 
ptr) = 2are r 2 @ (2-229 
= @ ro <.8@ 
and the APD or exceedance probability 
fe 
-Xro 
Prtr > ro) = 1 = Piro = © ro 2 @ €2-24) 
= 1 ro < @ 

How well the single parameter Rayleigh distribution fits 
observed atmospheric noise data can be seen in Fig. 4. The 
data points are exceedance probabilities plotted in decibels 
(dB) above the root-mean-squared (RMS) value of the received 
envelope. The data was measured on an ARN-2 receiver at a 
center frequency of 10 MHz at Boulder, Colorado [CRef. 13]. 
The high probability, low amplitude portions of the observed 
data approaches the slope of the Rayleigh distribution curve 


For « equal to 12. (a is defined in Eqn. 2-24) This value 
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Atmospheric Noise and Rayleigh Distribution 


25 


of ~« was chosen to fit the high probability portion of the 
data. That the data approaches a Rayleigh distribution is 
due to the effect of the many small overlapping impulses 
occurring at low amplitudes whose quadrature Gaussian 
components will have a Rayleigh envelope. However, at the 


low probability, high amplitude part of the curve the 


Rayleigh distribution predicts far fewer values than 
actually occur. This 1s the distinguishing factor of the 
envelope distribution of impulsive noise: that its 


distribution contains a higher probability of high amplitude 
terms than is predicted by the Rayleigh distribution. The 
second curve in Fig. 4 where «ais equal to .5 1s the 
Rayleigh distribution with the same power as the observed 
data. It can be seen that this curve under estimates the 
high amplitude and over estimates the low amplituae: Janie 
of the observed distribution. It should be noted that this 
type of plot emphasizes the low probability portion of 
distribution since that is where the deviation from the 
expected behavior occurs. 

In an attempt to correct the poor fit of the Rayleigh 
distribution at high amplitude levels, Likhter (CRef. 14] 
svaneued a combination of. two Rayleigh distributions: 

2 2 
“Ao -Bro 
Pr(r > rod = (l= eve + ce : (2-2ep 
This formula thas been shown to agree poorly with practical 


results (Ref. 15]. 
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Spaulding, Roubique and Crichlow CRef. 16] combined the 
Rayleigh distribution with a "power" Rayleigh distribution 


to obtain a distribution that fit very well for atmospheric 


noise over a wide range of receiver bandwidths. This APD is 
7a 
“Xo 
Pr(r > ro) = e rou « 6 
2i/s 
—- (Caro ) 
Pr(r > ro) = e Tome 6 C2=25) 


where a, 8 and 5s are determined from measured statistical 
parameters of the noise. 

Horner and Harwood (CRef. 17] proposed the two parameter 
log-normal distribution and found it gave a satisfactory fit 
to radio noise data in the VLF band. The log-normal was 
chosen since it has a more impulsive tail and fits the high 
amplitude, low probability data better than the Rayleigh 
distribution. The PDF for the log-normal distribution is 

2 2 
1 =CimCrn) .— a). /20 
De. = € (2-27) 
Lie 
ro(2n) 
Two examples of the APD of the log-normal distribution are 
shown in Fig. 5S where the curve with « equal to -.24 and ga? 
equal to .24 is the log-normal distribution with the same 
mean and mean square parameters as the curve with « equal to 


mo In Fig. 4. Comparing the two curves the log-normal 


distribution provides an intrinsically better fit to 
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impulsive noise data at both the high and low amplitude 
values. When the log-normal density has both parameters 
chosen to match the mean and mean square values of the 
observed data an excellent fit to this particular data is 
obtained. This 1S seen in the curve with «a equal to -1.4 
and g* equal to 1.4. 

Another model that has been successfully used to fit the 
observed APD of atmospheric radio noise is Hall’s general- 
ized "t" distribution CRef. 3]. This model is unique in 
that 1t 1s not of the filtered impulse type; however, in 
contrast to the strictly empirical models it does postulate 
a random process. 

One of the problems of Gaussian models of atmospheric 
noise is that the dynamic range of the model is less than 
the observed dynamic range ‘Ge the noise. To achieve a 
greater dynamic range starting with a Gaussian process, Hall 
proposed a model which considered the received noise to be 
of the form 

net) = acts (tb) C2=25)) 
where a(t) is a slowly varying stationary random process and 
s(t) 1S an independent, narrowband Gaussian process. An 
‘analytically tractable distribution was selected for a(t) 
which was chosen to give good agreement between the model 
and measured atmospheric noise data. The distribution for 


a(t) iss 


(ape, 


eae 
(m/2) 1 -(m/2a a ) 
3G = Eee 
m m+ 1 
o¢ I'(m/2) lal (2-22) 


where m anda are the two parameters used to fit the model 


to the data. The distribution of s(t) 1s Gaussian; 
2 2 
| -(s /201 ) 
“p(s) = —---———-_—C es ; 
e 1i/7e 
C2nq, ) (25a08 


Hall calculated the density function of n(t) to be 


sa 
r¢(8/2) Y 1 
ptn)d = err nea gee a, 
i 2 eevee 
lr ( (6-9) 72) Tt (n + y ) (2=31 
where 
ys 
Y=mMm (gq, /qa), 
8=me= il, 
and [Tf is the gamma function. For the case where a, = 0, 


EQn..ve-S! 1s the density function for Student’s "“t" 
distribution which is the basis for describing the density 
as a generalized "t" distribution. 
Hall then calculated the envelope distribution for n(t) 
based on the above assumptions and obtained 
g-1 r 
p Ge) a= Comal 7 
2 2 ((0+1)/2) 
Cr+ vee (2-32) 


with the phase uniformly distributed between @ and @2n. 
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For atmospheric radio noise, values of 8 in the range 2 
to 5 have been found to give excellent agreement with the 
data (Ref. LS) okos alee (Models with Rees vor values of the 
parameter 98 will be abbreviated to Hall2, Hall3, etc.) 

Fig. 6 shows the Hall3 fit to the same atmospheric noise 
data used in Figs. 3 and 4. An excellent fit to this data 
1s obtained. The Hall model thas some disadvantages. The 
higher moments only exist for orders greater than 6-1 and 
the parameters 8 and y must be determined for each 
interference scenario and are not easily related to the 


physical source of the interference. 


D. FILTERED IMPULSE NOISE MODELS 

The filtered impulse models differ fundamentally in 
concept from the empirical models described previously. 
This class of model will be explained in terms of the 
generalized interference scenario presented at the beginning 
of the chapter. The filtered impulse models have the 
advantage of being based on the underlying physics of the 
process) but suffer from the disadvantage of being 
analytically complex. Interestingly, some of the results 
obtained from the filtered impulse models have retroactively 
justified the expressions that were derived from the 
empirical models. 

To determine the APD of the envelope of a filtered 


impulse process, the joint characteristic function of the 
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inphase and quadrature components 1s determined. Since the 
noise 1S narrowband it will be circularly symmetrical and 
the 2-D Fourier transform will become a Hankel transform 
meet. 29,cl1l. This characteristic function 1s transformed 
to polar coordinates, then inverted to determine the APD. 
Starting with the complex noise process described in Eqn. 2- 
18, the joint characteristic function of the inphase ana 
Quadrature components 15 


j Cw. Ne Ca) + We Msg Ce) ) 
OCw, 9 Wa ) = ELe i (2-33) 


It 15 shown in Appendix B that the characteristic function 


of the envelope due to the Gaussian term is 


2 
=NoBCw. 72 
Da fwr) = e G2a=34) 


and the characteristic function of the envelope due to the 


impulsive term 1s 
w ah 


- 


sf pcan PoeGie aite Cy ea | idea C2=S5a0 





@ @ 
Or Cwy ) = e 


where pa) is the impulse amplitude density function, A, the 


pPemsison rate function, Me (bt), the inphase component of the 
filter impulse response and Jo 1s)6)6ltthe)06OCUlordinary)=—= 6—Bessel 
function of the first kind. The density function of the 


envelope is then the inverse Hankel transform given by 


ptr) = Pula (mi Oo CO) Oy CM. aw... C2=36>2 
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This equation for the PDF is difficult to evaluate; however, 
a number of researchers have examined special cases to 
obtain results. 

Furutsu and Ishida CRef. 22) derived an equation of the 
same form as Eqn. 2-36 with no Gaussian term. They obtained 
many approximate results by considering the equation for 
specific source density functions. If pCa) 1S exponential 
they found that the resulting envelope density was 
approximately Rayleigh for small amplitude values. If pla) 
was uniformly distributed, the log-normal density was found 
to be a good approximation over most of the range of 
interest. If a strong local source of interference was 
present, Furutsu and Ishida showed that a function involving 
the confluent hypergeometric Function was a good 
approximation. They also considered one special case were 
the unit counting function NCT) was modified to be a Poisson 
~ Poisson branching process and showed that for low 
amplitudes the resulting density function was approximated 
by the simple Poisson case. 

Giordano (Ref. 23] evaluated Eqn. 2-36 for various 
distributions of pla) which he determined from assumed 
propagation laws and spatial distributions of impulse 
sources. In one particular case he adopted the following 
assumptions: 


1. No Gaussian component 
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2. Uniform spatial distribution of impulse sources 
3. Inverse distance law for received field strength 
4. Arbitrary receiver envelope response 
and showed that the APD of the envelope function was approx- 
imated by 
Petr > re). = SS C237) 
2 2si72!) 


Cr. + K ) 


where 
t 


Ke | need at , 
@ 


T is the observation interval, h(t) the filter response, and 
cis ascenario dependent constant: This result is sig- 
nificant because it has the same form as the APD of the 
Hall2 model and physically justifies what had heretofore 
been an empirical oe. 

Middleton CRef. 24] included the Gaussian noise term in 
Eqn. 2-36 and expanded the characteristic function without 
taking the expectation. Then by ainverting the  char- 
acteristic function term by term, a canonical form of the 
envelope density function can be obtained. The exceedance 
probability was then shown to be 


Pror > Yaw 2 = 


m 
2 2a C-1) m ma ma 
“ro C1 - ro DCm AL POL + LF C1 - 3 23 1) ) 
m=O m'! 2 2 


(2-38) 
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where ,F, is a confluent hypergeometric function, and a and 
A, are two parameters that are determined by the source dis- 


tributions and the propagation law. 


E. WALL MODEL EXTENSIONS 

For later work in this dissertation a model “for 
atmospheric noise statistics is required. This section will 
further describe the Hall model as applied to HF atmospheric 
noise and will derive some extensions to it that will be 
used to digitally simulate atmospheric radio noise. 

The Hall model for the envelope statistics of radio 
noise 1s 

= 1 r 
ptr) = (6 - 1) vy 
2 2meee 1) 72) 

Cr +y ) (2>=3S9 
and was chosen for a number of reasons. The most important 
are: the Hall model is analytically tractable and shows a 
good fit to the data. It can be obtained for certain values 
of the parameter 6, starting from the physically based 
filtered impulse models. This was shown by both Giordano 
CRef. 23] and Middleton CRef. 24]. Additionally, Schonhoff 
CRef. 18] has related the Hall models to CCIR Report 332 
CRef. 25], the standard reference of first order statistics 
of atmospheric noise. In the CCIR report the APD curves are 
shown as a function of Vo. The parameter, Vo, 1S a commonly 


used measure of the impulsiveness of atmospheric radio noise 


and is defined as 
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Bl Mf 
a 
Vo = eG@log 


Elid Ce=40) 
where r is the amplitude of the received envelope. For 
Quadrature Gaussian noise Vp = 1.05 dB and as the noise 


envelope becomes more impulsive, Vy will increase. 
One problem in relating the Hall model to Vn 185 that the 


moments of the Hall model only exist for orders greater than 


a |. This means that for the Halle and Hall3 models, Vp i185 
undefined unless the model 1s modified. Two methods have 
been proposed to do this. Hall’s method (Ref. em adds 


another term to the envelope density function and multiplies 
by a negative exponential term to reduce the tail of the 
density function and give finite moments. The advantage of 
this method is that the density function exists from @ to 
mat initcy ; however, the resulting density function 15 
complicated. The method proposed by Schonhoff (CRef. esi a: 
truncates the Hall2 and Hall3 distributions above a set 
level, Te, and renormalizes the density function by picking 


kK such that 


Te 
ke pervar = 1. 2-4) 
0 
The level, Te, 15 chosen to set a desired Vp for the 
Bes tie 2 DUt1ON. The constant k will be a function of OD where 
er il/e 
Cie Ole /yV) ) (2-42) 


Sy, 


The existence of the first and second moments for these two 
distributions 1s then assured. Using this’ procedure 
Schonhoff generated a family of “ieee ueneene parameterized 
by 6 that can be used to represent a wide range of measured 
data. By matching the Vp ratio of the data to the dis- 
tribution, a close fit to the APD is found. 

An extension to the basic Hall model to be used ina 
later chapter 1s the density function of a constant 
amplitude sinusoid plus the assumed noise. This density 1s 
analogous to the Rician density in Gaussian noise theory. 
It will be used in simulating the performance of energy 
detection receivers'7 for. the signal plus noise case. It was 
derived in (CRef. 26] for the Hall2 model and it is derived 
for the Hall3 and HallS models in Appendix C. For the Halls 
model the density function 1s 

2 2 2 2 
(r + AWt yy ) e2y 7D 
or) Fe eee 
4 22 Ze eve 2 4 2 
(r - GA r + @2A Y 2+ @y Fetes) yu oon 
where A 1S the signal amplitude, y 1s a Hall3 noise 


parameter and Dis’ defined in Eqn. 2-42. Likewise the 


density function for the HallS model is 


a 2 2 

Cy r(b + 2a ) 

Te (2-44) 

2 e S/e2 

Ca - b) 

where 
2 2 2 
a= r +A + ¥y 
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and 

Bo > -2Ar. 
Fig. 7 is an example of the APD for a sinusoid plus Halls 
noise with the RMS value of the noise set to 1 and the 
Signal to noise ratio at -26, -6 and 14 dB. For the -e6 dB 
SNR case the APD is very similar to the Hall3 noise only 
density shown in Fig. 6. 

Commonly used functions relating to the Hall models are 
shown in Table I. These are the envelope density function, 
the lst and 2nd moments of the envelope, the inverse of the 
CDF, the quadrature marginal density function and the Hall 
noise plus random sinusoid density. The parameter D in 
Table I is defined by Eqn. 2-42 and the function €C€:J is the 
complete elliptic integral of the second kind. The quad- 
rature components can be derived by transforming the 
envelope density to rectangular co-ordinates and integrating 
with respect to one component. The inverse of the envelope 
CDF is also included since it provides an easy method to 
generate random deviates. 

Also useful for simulation purposes are the values of 
the Hall parameters, 9 andy, for different values of Vp 
normalized to a unit root mean square value. The truncation 
point Tp was found for a desired Vy using an iterative tech- 


nique. These parameters are shown in Table II. 
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Ve Hall No. Elm toa cine) 2D Tp Y D 
2.1 Ss, 1 s7oo = 1.@1 = 

3 4 1 . 708 = - 708 = 

4 3 4 =o3t 11.66 - 420 277s 
=) 3 1 - 62 25.90 - 364 71.05 
6 3 1 -201 G7 eae soe. 2 ies! 
8 2 1 ols se ye - 985 137.6 
1@ a 1 oe 17.08 - 058 291.0 
2 e 1 sem t 24.25 - 041 S97 
14 2 1 mh, 33.81 "oes 1141.9 

TABLE II 


F. SUMMARY 

Both empirical and analytical filtered impulse models 
have been successful in explaining the first order envelope 
statistics of bandpass radio noise processes. However, in 
obtaining these results the assumption has to be made that 
the impulse arrival times due to the noise process form a 
Poisson point process. In general this will not be true for 
man-made impulsive noise and for some types of atmospheric 
noise (Ref. 27). In most cases of man-made noise the 
underlying deterministic mechanism of the device causing the 
noise will modulate or regulate the impulse distribution. 


The Hall model for atmospheric noise has been examined in 
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detail and some useful extensions to it have been derived 
that will allow noise and signal plus noise environments to 


be simulated. 
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III. GAP NOISE MODEL 


A. INTRODUCTION 

In this chapter a specific source of man-made radio 
noise will be modeled and specified in terms of the counting 
process, N(t), and the probability density function (PDF) of 
the impulses, ptfa)d). These specifications were introduced in 
the previous chapter to define the impulsive noise model. 
The noise source that will be examined 1S gap type 
discharges and 1S commonly found on electric power dis- 
tribution and transmission lines. Gap noise 1s one of the 
major types of interference from power lines and is 
frequently observed as the primary interference to com- 


munication systems operating in the high frequency (HF) 


radio band. 


B. GAP NOISE MECHANISM 

At least two mechanisms have been found by which a gap 
discharge process can occur on a power line. The resistance 
in the line insulators can be degraded allowing current to 
Flow through the insulator base and creating a potential 
Gradient across any gaps or defects in the insulator 
mounting hardware. A second way in which a potential can be 
created across anair gap 1s by an electro-static coupling 
of the line potential to isolated hardware on the pole. In 


both cases the potential across the gap is discharged by the 


a 4 


voltage breakdown of the gap and the resulting rapid current 
flow or spark. This process generates a radio frequency 
CRF} noise impulse with spectral components extending into 
the hundreds of MHz. During a Single discharge, the 
potential across the gap 1S temporarily diminished. WHow- 
ever, while the fundamental 60-H2z waveform is still greater 
than a threshold voltage, the process can occur again. The 
spark will discharge across the gap . repeatedly until the 
alternating current waveform drops below the breakdown 
threshold potential. CRef. 28:p. 78,29] 

The spark discharge and recharging of the gap potential 
indicates that this type of process 1S regenerative and can 
be modeled as a renewal process where the renewal points are 
associated with the sparks. One feature of gap noise that 
complicates the modeling is the 1290-Hz on-off-on modulation 
imposed on the renewal points by the alternating current 
CAC) waveform. One way to account for this effect is to 
consider the turn on times as another renewal process driv- 
ing the spark discharge process. This type of model 1s 
known as a branching renewal process (Ref. 30]. 

The noise processes that were used to develop this model 
were short term stationary in that the statistics of the 
process did not change over the observation interval. Not 
all sources of power noise are stationary even over a short 
observation interval. Physical effects such as wind, solar 


heating and varying line loads can act to make certain types 
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of power line noise highly variable from observation to 
observation. Although these types of nol1se were observed 
during data collection, the data used to develop the model 
was taken from noise sources that were stationary for the 


length of the ten minute data records. 


C. OBSERVED TIME-DOMAIN CHARACTERISTICS 

The gap noise sources for this study were observed on 
utility distribution lines in the vicinity of the Naval 
Postgraduate School and were chosen to illustrate parameters 
of the noise model. Fig. 8 1s a typical time-domain observ- 
ation of agap discharge process observed at 3 MH2 = and 
envelope demodulated with a 1@0-kHz Gaussian bandpass filter. 
The important characteristics to note are: 

1. the process has an on-off-on modulation at a 120- 
Hz rate related to 6@0-Hz waveform of the power 
line, and 

2. the pulse groups that result from the modulation 
have a probabilistic number of impulses occurring 
in each group and a probabilistic interarrival 
time between pulses 1n a group. 

In the first pulse group, 9 pulses occur with varying 
amplitudes and interpulse arrival times. In the second 
pulse group starting approximately 8.33 msec later, 8 pulses 
occur, agaln with varying amplitudes and interarrival times. 
One group of pulses is associated with either the positive 


or negative polarity of the line voltage waveform and the 


other group with the opposite polarity. Identification of 
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Figure 8. Envelope of Demodulated Gap Noise 
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the polarity 1s impossible without physically locating the 
source. 

The average amplitude of the impulses will be one of the 
parameters for the model. An estimate of this parameter can 
be made from the above presentation. It is important to 
note that for this highly impulsive type of noise that the 
observed amplitude is a function of the shape and width of 
the effective bandpass filter and the detector char- 
acteristics. 

At this point it should be noted that the model can be 
approached at two levels of complexity. In the simple 
version the data from both the positive and negative phases 
1s considered as a whole and averaged over the two phases to 
determine a set of average parameters. In the more complex, 
but more accurate model, two-sets of model parameters are 
determined - one for each phase. For the remainder of this 
work the more complex case will be considered. The results 
for the simple case are presented in Moose and O’Dwyer (CRef. 
31] and can easily be determined from the complex case by 
setting the parameters equal for both phases and simplifying 
the resulting equations. 

In order to better characterize the interpulse arrival 
times, which will be used to determine two model parameters, 
the rising raster capability described in the § instru- 
mentation section was used to generate the display of a 


different gap noise source which is shown in Fig. 9. In 
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Figure 9. Multiple Scans of Gap Noise Interference 
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this picture the amplitude of each individual record was 
reduced and thresholded so that all that essentially remains 
is the time of impulse arrival data for 28 time records. 
The time base of the display was intentionally synchronized 
to the power line to facilitate taking data and this 
accounts for the regularity of the pulse groups’ from 
observation to observation. In this view the interarrival 
times for S56 pulse groups (28 groups of one polarity and 28 
of the opposite polarity) can be determined along with the 
number of pulses in each of the S56 successive pulse groups. 
The average number of pulses per group considered separately 
for each phase will also be used as estimates of two model 
parameters. 

Fig. 10 is a histogram of the distribution of the 
interarrival times between the observed impulses for the 
Same noise process shown in Fig. 9. To use this data to 
generate parameters for the noise model, the histogram will 
be approximated by a continuous density function. The gamma 
density function CRef. 32e:p. 18] 


r-1 -CAt) 
p(t) = ACA) e /(r-1)! (3-1) 


was chosen because it showed a close fit to the data, and 
its characteristic function, which will be used in later 
derivations, was particularly simple. The characteristic 
Function of the gamma density is 


Se 
OC 5005 —— C1 = wr : (3-2) 
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The r andA parameters are simply interpreted in terms 
of an underlying Po1lsson process as the time tomiee misin 
point of a Poisson process of intensity A. If wp is the 
random variable assigned to the observed interarrival times 
then the unbiased estimate of the mean of p 15 

= N 
est(yp) = €1/N) £ py (32a) 
i=1 
and the unbiased estimate of the variance 15s 
e 


est(o, ) = ¢ 
1 


2 _ 2 
Hi - N estCp) }/¢(N-1). (3-4) 
1 


m2 


The unknown parameters of the desired gamma density function 
are estimated using the method of moments [CRef. 33:p. 250] 
using the following relations; 


es — 2 a 2 
r= est(p) /estla ) and A = est(pd/est(aq, ). CS>or, 


Using this approach, the gamma density function used to 
approximate the interarrival time histogram is also plotted 
on Fig: 10. Figs. 11 and 12 are based on the same data as 
Fig. 19. Wowever, they are separated according to different 
phases. It can be seen from the r and 2A parameters in Table 
III that there 1s a significant difference between phases of 
the same gap noise process shown in the histogram even 
though the mean value of each data set was almost the same. 
The data in Fig. 11 was more clustered than the Fig. 12 


data. This results ina higher r and lower 2X parameter for 
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the 


gap 


the 


Case 


Case 


Case 


Case 


Case 


A F1 


same mean value. PaAgs. 1G ands!14 are additional sets of 
noise data that were observed during the research. 
Table III summarizes the rand A parameters fitted to 


data described above. 


r x 
il SS.c 37,909 
1 Phase A 63.4 71,400 
1 Phase B 23.3 25,9900 
2 44.3 74,600 
S 14.1 19,2900 
TABEEVITI 


nal model parameter will be called Tp and 1s a constant 


delay or offset of every second pulse group that exists with 


respect to the fundamental voltage waveform. 


data 


In summary the inputs to the model based on time-domain 
are; 


1. An estimate of the average amplitude of all 
observed pulses. 


2. An estimate of the average number of pulses per 
group for each phase. 


3. An estimate of the mean of the interarrival time 
of the pulses for each phase. 


4. An estimate of the variance of the interarrival 
time between pulses for each phase. 


5S. An estimate of the delay in the start time of the 


pulse groups of one phase with respect to the 
fundamental voltage waveform. 
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The average amplitude of pulses from different gap noise 
sources is highly variable, ranging from the instrumentation 
noise floor to higher than any observable signal in the HF 
band. The average number of pulses per group has been 
observed to vary from 1 to greater than 2@. The mean of the 
interarrival times ranges from el to 1 ms and the 
coefficient of variation (standard deviation divided by the 
mean) of the interarrival times has ranged from approx- 
imately .1 to .3 for the gap noise processes we have 


observed. 


D. MODEL DEVELOPMENT 
The principal goal of this chapter is to describe a 
noise model that, with suitable choice of parameters, pro- 
vides an adequate simulation of the actual physical noise 
mechanism of gap noise. In this section the time-domain 
data and parameters derived in the previous section will be 
integrated into the filtered impulse model developed in 
Chapter II. Recalling from Chapter II, a filtered impulse 
noise process 1S given in complex envelope notation by 
NC) = 40% 
n(t) = g(t) + LF ae RG ee ee (3-67 
i1=1 
where pCa) 1s the amplitude distribution of the impulses, 


N¢Ct) 15 a unit counting process whose statistics determine 


the impulse arrival times, 8 is a random variable uniformly 
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distributed between @ and en AG « co 1s the impulse 
response of the receiver filter. 

In order to physically justify a filtered impulse model, 
the impulse duration must be small compared to the inverse 
bandwidth of the receiver filter. This condition 1s easily 
met in the case of gap discharges. Laboratory analysis of 
temporal characteristics of gap discharges for various 
geometries shows that the impulse durations range from 19 to 
190’s of nano-seconds [CRef. 34]. Therefore, tOn wat lt cer 
bandwidths up to 1 MHz the output noise process will only be 
a function of the incident time of the impulse and the 
filter response, not the waveform of the impulse. 

The specification of the statistics of N(t)d, the count- 
ing process that drives the model is then crucial to 
obtaining an accurate representation of the physical noise 
process. Based on the electrical characteristics of a gap 
discharge and the observed time-domain characteristics, the 
Following assumptions are made; 

1. A primary series of event times separated by an 
interval Te/2 exists. To 15 the fundamental 
period of the power line voltage waveform. The 
distribution of the time to the first primary 
event is uniform over @ to To. 

2. A subsidiary process commences at the primary 
evene-cimes 1.6. lo/2 5 lo @eSsle/e2 ote The 
subsidiary process 1S a renewal process that 
continues for Ni or N2@ renewals alternating 


between Ni and N2@ renewals at successive primary 
Dolnmcs. 


SPS, 


3. The interarrival times to the first and subsequent 
points of the subsidiary renewal process for each 
Phase are all independent and identically 
distributed CIID) within a pulse. Each phase has 
its own set of parameters. 


4. The weighting or amplitude distribution of all 
impulses in the subsidiary process 1S a constant. 


5S. A constant offset or delay designated by Tp may 
exist at every second primary renewal point. 


The first assumption is supported by the periodicity in 
the data that is related to the fundamental frequency of the 
power line. The large majority of observed gap noise sources 
had pulse groups on both the positive and negative phases of 
the fundamental waveform. This accounts for the T./e per- 
LOdici£y. 

The second assumption addresses the differences between 
the sparking phenomenon on the positive and negative phases 
of the fundamental. It 16 frequently observed that the 
ECNeos] is different than ECNwea J] where Neos and Nwca are 
the random number of discharges in the pulse groups. This 
can be due to two effects; an asymmetrical gap geometry and 
the fundamental physical difference between the sparking 
mechanisms for positive to ground and negative to ground 
sparks. Two deterministic integer constants; Ni and Ne 
nearest to ECNeosJ and ELNvee 2 respectively are assumed as 
the approximation to the variable number of sparks per half 


cycle of the fundamental waveform. 
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The third assumption is based on the fact that the 
“inception of gap discharges in natural air and the 
development of electron avalanche are fundamentally 
probabilistic processes that depend on atmospheric pressure, 
humidity, presence of natural ions, electrode surface and so 
on" CRef. 34]. In view of the above statement and 
considering the empirical data, the justification for devel- 
oping the gap discharge as a probabilistic process is well 
founded. The assumption of independence from discharge to 
discharge 1S not as well justified. Effects such as 
electrode heating after the initial discharge ina cycle 
could act to make the average interarrival time vary from 
pulse to pulse within a pulse group. However, to develop a 
tractable model the assumption of independence between 
arrival times 1s made. The gamma density function deter- 
mined by the estimate of the mean and variance of the 
interarrival times for each phase is used to define the PDF 
of the interarrival times. 

The fourth assumption concerns the amplitude of the 
impulses. In most cases the amplitude of the impulses 
within a pulse was nearly constant. In some cases there was 
a difference in the average amplitude from negative phase to 
positive phase, however to keep the model tractable the 
amplitude of all impulses 1S assumed constant. 

The fifth assumption incorporates an effect that will 


explain one of the features observed in the power spectral 
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density (PSD) of the observed process. The physical 
rationale behind this assumption 1s that there may exist a 
difference i1n the threshold between the positive and 
negative phases, thus consistently delaying or offsetting 
the start of the renewal process on one of the phases. 

The above assumptions place this model in a class of 
processes known as branching processes with the main process 
being a degenerate renewal process with an interarrival PDF 
of 

p(t) = &¢(t-T, /2) (327) 
The subsidiary process is a renewal process of Nl or Ne 
points. Consequently, the impulse source for the gap noise 
model is completely specified by the following parameters: 


a - Amplitude of the impulses, 


N1,N2 - Number of impulses in the subsidiary process, 
HWi,He - Mean value of impulse interarrival time, 
Vardyp,),Var(pe) - Variance of interarrival times, 

To - Of Fset on one phase. 


To complete an interference scenario NoBerr /2 and H(t) must 
also be specified. The expression for the branching renewal 
process is then 
MCT) Ne Oe 
nt) = ¢ (tees , Aen @ nee ce a lee (3333 
m=1 n=1 


Fig. 15a,b shows sample realizations of the envelope, 


ECt), of the noise process defined above. 
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Figure 15. Envelope and Counting Function of nt) 


The Gaussian term was set to @ and H(t) was approximated as 
a decaying exponential. Note that in comparison to an 
actual noise process illustrated in Fig. 8 the amplitude of 
the impulses is constant with a fixed number of impulses on 
the positive and negative polarity of the line voltage 
waveform. Fig. 15c is the primary counting process, M(T), 
driving the subsidiary process shown in Fig. 15b. Fig. iSb 
includes a non-zero Tp and its relationship to MCT) is shown 
moe ig. iSc. 

eG. 16 shows a realization of the envelope of this 
process for Nil and Ne = 3, T = 30 ms and gamma density 
function parameters rl and re = 32 and Al and A2 = 37,9000. 


Random interarrival times with a gamma distribution were 
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approximated Dy summing exponentially distributed random 
variables [Ref. 35]. The impulse amplitudes were normalized 
to one anda background level of Quadrature Gaussian noise 
at NoB/2 = .9@0@1 was added which simulates either receiver 


noise or high density, low amplitude impulsive noise. 


E. COMPARISON OF CALCULATED AND OBSERVED DATA 

One partial description of a noise process 18S the 
amplitude probability distribution (APD) of the’ envelope. 
This description has been used extensively in the analysis 
of noise and in deriving optimum receivers. For the N(t) 
specified in the previous section this calculation would be 
errr cult. A second noise process descriptor 1s the 
spectrum of the noise process. This descriptor is suited to 
our model where the structure of the noise 1S contained in 
the counting process that drives the impulse generation. In 
addition to receiver noise performance evaluation, spectral 
analysis can also be used for noise source identification 
and isolation (Ref. 36]. 

In order to determine the spectrum of the envelope 
squared of the branching renewal process defined by Eqn. 
3-8, the spectrum of the envelope squared of a non-branching 
renewal process is determined in Appendix D. The results 
obtained in Appendix D are then used to derive the spectrum 
of the the branching renewal process which 1s postulated to 


model gap noise. Starting with the branching renewal 


SS 


process given in Eqn. 3-8, and assuming g(t) is equal to 
zero, it 1S Shown in Appendix E that the average of the 
Bartlett estimate of the PSD is: 
= S 
Scw,T) = |He(w)| /T © Term (Cn) (3-9) 
n=1 
where n 1S an index to the 5 Terms defined below. These 


five Terms involved in the summation in Eqn. 3-9 are defined 


as follows: 


Term 1 - The inter-pulse group summations between pulses in 
pulse groups with Nil pulses 

Term 1 = M2 2(w,Nl1) (3-9a) 
Term 2 - The inter-pulse group summations between pulses in 
pulse groups with Ne pulses 

Term 2 = M2 Ss (w,Ne2) (3=96> 
Term 3 - The intra-pulse group summations between pulses in 
groups with Nl pulses 


cos(M2 2nw/w.e) - 1 


Term 3 = € rrr" eM FF 2Cw,N1) (3-S9e 
cos(2nw/w.) - 1 
Term 4 - The intra-pulse group summations between pulses in 


groups with Ne pulses 


cos(M2e 2Cnw/we) - 1 


Ten 4) Aten eoeqre_—X—X—X——a-__—————( r:iéCOM! Nw, NA) (3- 9a 
cos(Cenw/wo) - 1 
Term 5 - The intra-pulse group summations between pulses in 


groups with Nil pulses per group to pulses in groups with Ne 


pulses and vice versa 
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cos(Me enw/we) - 1 


Term 3S = @cos((wlT./2-lh)) (errno Ft By oN1 De 5 N2D 
cos(2nw/w,) - 1 

(3-9e) 

where @, and @2 are the characteristic functions associated 

with the interarrival times for each phase. The parameters 


N1 and Ne are the number of renewals associated with each 
phase and @&, 2, and w are defined below: 


2 N 
N - NI@¢cjw) | OCjw) (1 - Olyw) dD 
S(Q,N) = ———————— _ -- Ss 2? Ref ——————————__———_ }, 
2 2 
[1 - @¢Gw)| (1 - B(5w)) (3-10) 


N+1 2 
[SC jw) - BOCjw) 
2(0,N) = ———__________——__ ,, 
2 
}1 - @Cjw)| (3-11) 


and 


Ni+t+1 Ne+1 

(O@, (jw) - BO, (jw) > (B82 C-jw) - BO. (-jw) ) 

wtO, ,N1,02,N2) = : 
C1 - @, (gwd) ¢€1 - Bea C-gw)) 





C3129 

As a test of the model the PSD of the gap noise process 
shown in Fig. 9 was determined using a spectrum analyzer and 
also analytically determined using only the model parameters 
from the time-domain data and the methodology of the 
previous section. Fig. 17 1S a computer plot of the 
analytic estimate of the PSD which is compared with the 
observed PSD computed on a Wavetek UAS@@A spectrum analyzer 


shown in Fig. 18. 
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Figure 18. Observed Estimate of Power Spectral Density (PSD) 
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Figure 17. Predicted Power Spectral Density (PSD) 


Comparing the two PSD’s, it 1S seen that the analytical 
expression correctly predicted all of the significant 
features seen in the observed PSD. 

The components of the spectrum due to the different 
terms is of interest and offers some insight into the 
origins of the features of the observed spectrum. The dom- 
inant feature of the observed spectrum are the periodic 120 
Hz spectral lines. At 800 Hz the spectral lines change, and 
appear as odd harmonics of 60 Hz, however, still remaining 
at 120 Hz intervals. Terms 3 - 5 are responsible for the 
these spectral line features. Each of these terms is a comb 
function, determined by the fundamental frequency = and 


observation interval, multiplied by an envelope determined 
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Figure el. 


by the interarrival time distribution. Figs. 19 and 2@ are 
Terms 3 and 4 -respectively and clearly show the phase to 
phase difference. The more peaked behavior of Fig. 19. at 
1100 H2 is due directly to the more clustered interarrival 
time behavior on one phase. This effect was observed in the 
time domain in Fig. ile. When the terms shown in Figs. 19 
and 2@ are added together, they form 60 Hz harmonic spectral 
lines and when Term 5S, shown in Fig. 21, is added, positive 
reinforcement occurs at 120 Hz and negative at 60 Hz 
harmonics. This produces the observed 120 Hz harmonic 
spectral lines. When the cancellation is not complete due 


to differences in the characteristic function from phase to 
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Figure 19. PSD Term due to Phase A Intra-Pulse Group 
Correlations j 
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Figure 20. PSD Term due to Phase B Intra-Pulse Group 
Correlations 
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phase, small 60 HZ components can be present. This effect 
1S seen in both the observed and predicted spectra. 

The observed spectra changeover from 120 Hz harmonics to. 
odd 6@ Hz harmonics at 8@@ Hz 1S also correctly predicted by 
the model. This effect 1S due to the factor, To, appearing 
in Term 3S which causes the odd 6@ H2 harmonics to be 
positively reinforced and the 120 H2 harmonics to be 
Canceled out above 8@@ Hz. The final effect shown in the 
model and seen in the observed spectra is the presence of 
continuous spectral components. The continuous spectra due 
to Terms 1 and 21S shown in Figs. 22 and 23 and =«also 
differs from phase to phase. The continuous spectra peaks 
markedly at 1100 Hz and this effect 1S apparent in the 
observed spectra as aA rise in the note floor. The come 
tinuous spectral term arises when Term 1 1s added to Term 3 
and Term eto Term 4. The peak value at 1100 Hz in Terms 1 
and 2 is greater than the negative values in Terms 3 and 4 
at the same frequency and when added together causes a 
continuous rise in the spectral floor. The reason for this 
1s that the 1 to i+1l and i-l interarrival times within a 
pulse group are more correlated than the corresponding 
interarrival times between different pulse groups. The more 
correlated arrivals appear in Terms 1 and 2 and cause the 


spectral peak to be larger. 
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Figure 22. PSD Term due to Phase A Inter-Pulse Group 
Correlations 


MAGNITUDE 


900.6 600.0 


400.0 





SOSH o0.0 750 ,Gmn000.0 lec. OmIS00.0) 7500) 2000.0) 2250-0 
PREQUGNEY =- HZ 


Figure 23. PSD Term due to Phase B Inter-Pulse Group 


Correlations 


Zea 


F. SUMMARY 

A nine parameter model for a single source of gap noise 
interference was He venceee and shown to accurately predict 
the PSD of narrow bandwidth envelope demodulated gap noise 
interference. If less accuracy 1S desired, phase to phase 
differences can be neglected and the resulting model has 
only five parameters. The model was specified in terms of 
an impulse driving function, and techniques to estimate the 
model parameters from time domain data were described. 
Since the model specifies a complex driving impulse function 
it can be used aS a Simulated interference source for 


arbitrary filters and receivers. 
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IV. CORONA NOISE MODEL 


RE 


A. INTRODUCTION 

In addition to gap noise which was discussed in Chapter 
III, another major source of power line radio noise 1s 
corona noise. Corona is different in many respects from gap 
noise, and the model used to describe it will take a 
different approach. In contrast to gap noise, which is 
observed on both electric power transmission and dis- 
tribution lines, corona noise is observed only on power 
transmission lines and generally has a fundamental frequency 
of 180 Hz corresponding to the three phases of a 60 Hz 
alternating current (AC) system. Corona noise is caused by 
a partial breakdown of air surrounding a conductor which is 
at a high potential. Consequently, the impulsive structure 
of corona noise 1S not as well defined as the highly 
impulsive structure of gap noise. Corona noise from AC 
sources appears somewhat like amplitude modulated Gaussian 


noise. 


B. CORONA NOISE THEORY 

When an increasing potential difference 1s applied 
between two electrodes, a breakdown voltage 1s reached that 
1s characterized by the transition of air from a poor 
electrical conductor to a good conductor. If the field in 


the gap between the two electrodes exceeds the electric 


Uz, 


strength of air, which 1s about 30 kV/cm for gaps over 1 cm 
in length, prior to the spark Dreakdown voltage, then corona 
will axe Sharp points and small radius wires are avoided 
whenever possible on power lines, and corona generally does 
not become a significant source of radio interference until 
voltages exceed about 7@ kV (Ref. e28:p. 84]. 

As the voltage on an AC line is increased, corona will 
appear and be associated with either the positive, negative 
or possibly both half-cycles of the fundamental waveform. 
For an aged, clean transmission line it appears that the 
threshold voltage for negative corona is less than that of 
positive corona. If the transmission line has any unusual 
characteristics, the threshold voltages of positive and 
negative corona become difficult to categorize (CRef. 37]. 
In particular, it has been observed that positive streamer 
corona, which 18S more disruptive of communications’ fre- 
quencies than negative corona, predominates on newly 
installed lines and during periods of precipitation. For 
modeling purposes, the important characteristics of corona 
noise are not a precise description of the physical 
mechanisms but the relationship between the observed inter- 
ference waveform and the line voltage. 

The procedure most commonly used for the prediction of 
interference due to corona noise gives the noise power as a 
function of the characteristics of the interference 


scenario. The random properties of the corona noise are 
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averaged out CRef. 38]. These characteristics would 
commonly include distance from the line, the line voltage 
and a weather correction term. Each of these terms 15 
experimentally verified and related to some physical process 
in the generation of corona noise. While this procedure 
Gives an estimate of the noise power in a specified 
bandwidth, 1t does not incorporate the probabilistic nature 
or time-domain periodicity observed in the noise. 

The approach taken in this chapter will continue with 
the filtered impulse theme of the dissertation. For corona 
noise, the counting process, N(T), will be modeled as a 
periodically modulated stream of impulses. This approach 
will have the advantage of being accurately able to simulate 
the time- and frequency-domain behavior of corona noise. 
However, as with the gap noise model it will suffer from the 
disadvantage that for each scenario, a set of parameters 


must be experimentally determined. 


C. TIME DOMAIN CHARACTERISTICS 

The time-domain data for corona noise was obtained using 
the instrumented van described in Appendix A. The antenna 
was a S3-meter whip mounted on the van (which was’) parked 
underneath a 3-phase SO@-kV transmission line). The weather 
conditions at the time of measurement were clear and sunny. 
Records of the corona noise were obtained using the envelope 


detected output of the HP-141T described in Appendix A and 


TD) 


recording the data on a 2@-kHz bandwidth, 85-dB dynamic 
range digital audio system. The recording process allowed 
subsequent analysis in the laboratory. 

A 37-ms portion of a data record is shown in Fig. 24. 
The data was taken at a 30@0-kHz center frequency and 190-kHz 
intermediate frequency (IF). This particular picture was 
obtained by playing back the tape, re-digitizing the data 
and displaying it with a signal-analysis software package. 
One of the important features of the time-domain record in 
Fig. 24 1s the periodicity. The fundamental frequency of 
18@ Hz can be determined from a periodic rise in the mean 
value of the noise and from a periodic increase in high- 
amplitude impulses. The fundamental noise frequency of 18@ 
Hz on a three phase line indicates that the dominate noise 
1s generated on one half-cycle of the fundamental 60-Hz 
waveform on each of the three phases of the line. This 1s 
consistent with low-level negative corona predicted for aged 
lines 1n clear weather. 

The periodogram spectrum of the data record described 
above is shown in Fig. 25. The time-domain observations of 
the 18@-Hz periodicity are confirmed and a 360-Hz harmonic 
of the fundamental is also present. At higher frequencies, 
no harmonics are observed and the spectrum becomes white. 
This 1S in distinet contrast to the gap noise described in 
Chapter III where the spectrum of gap noise has harmonics 


well into the kHz range. The periodicity can also be seen 
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in Fig. 26 (which 1s the same data recorded on the 3-Axis 
display). The contribution of each phase of the power line 
to the total corona noise process can be clearly seen in the 
lower picture in Fig. 26. 

Another important feature of the data for modeling 
purposes is that the receiver response to ani individual 
impulse overlaps with the response of preceding impulses. 
Although the noise is still impulsive, the simplifying 
assumption of non-overlapping pulses 1s not valid. Leis 
seen in Fig. 24 that the amplitude of the impulses is random 
and there 1s no fine-grain interarrival time structure other 
than the periodicity discussed above. Both of these 
observations are in direct contrast to gap noise where the 
amplitude from impulse to impulse is almost constant, and 
there is a definite impulse interarrival time structure. 
These observations provide a method to distinguish between 
gap noise and corona noise. 

An interesting aspect of noise from transmission lines 
was noted while taking measurements for the corona model. 
It was difficult to obtain a recording of power-line noise 


in the vicinity of high voltage transmission lines due only 


to corona noise. In many cases, a gap noise component was 
superimposed on the corona noise. This can be seen by 
comparing Figs. 26 and 2/7. Fig. 27 1s a view of both gap 


and corona noise. The corona noise amplitude peaks at the 2 


yV level, and the amplitude of the gap-nolise impulses 15s 
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Corona Noise as Shown on 3-Axis Display 
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Gap and Corona Noise 


Figure 2/7. 


Ss 


about 3.5 wv. The two lower pictures are of identical data 
with the threshold varied to emphasize each noise type. The 
lower picture has the threshold set at about 3 wypV 
eliminating the 18@-Hz corona noise and showing the _ gap 
noise process occurring at a 12@0-Hz rate. This corresponds 
to a gap noise source sparking on the positive and negative 
half-cycle of a single phase of a 3 phase line. When the 
threshold is lowered to about 1 pV, the 180@-Hz corona noise 
dominates as seen in the middle picture of Fig. 27. This 
phenomenon 1s due to corona discharge on the negative half- 
cycle of each phase of a 3-phase system. Without the abil- 
ity to examine the time- and frequency-domain behavior of 
the noise using a scanning analyzer and 3-Axis display, it 
would be very difficult to differentiate between gap = and 
corona noise and make accurate measurements of either. 

Fig. 28 illustrates corona noise with different phase- 
to-phase characteristics. In this example the noise power 
due to the corona 1S unequal from phase to phase. The noise 
occurring on the dominant phase 1S approximately twice the 


peak voltage of the phase with the lowest amplitude. 


D. CORONA NOISE MODEL 

Based on the above observations of the demodulated 
envelope of corona noise, a specific filtered impulse model 
based on Eqn. 2-18 will be developed. This model 1s similar 


to models used in synthetic hydrology to account’ for 
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seasonal variations in streamflow data (CRef. 39]. Assuming 
no Gaussian noise component, a lowpass complex filtered 
noise process can be written 
NCt) - 59,4 
n(t) = Cae Ou Ce = tags C4-1) 
i=1 
Since corona noise is a high density process, the envelope 
cannot be simply defined as in Eqn. 2-22 since the response 
of the filter to adjacent pulses will overlap. To “Gites 
cumvent this difficulty it will be assumed that the envelope 
of n(t) can be modeled by 
Nct) 
i=1 
where both e(t) and g(t) are real, positive and g(t) depends 
on h(t) and the process intensity. : 
The assumptions to define the model are as follows; 

1. The rate parameter for the counting function N(t) 
1s periodically modulated at frequency wWo and will 
be approximated by a finite number of terms ina 
Fourier series. 

2. The amplitude distribution of successive pulses 
Given by the a ‘Ss are independent and identically 
distributed with density function pla). 

3. The start time of the periodic modulation is 
uniformly distributed over the period T. where To 


= 1/Wo ° 


Assumption one is motivated by the periodicity observed in 


the data. In most cases the fundamental period, wWo will be 
18@ Hz. If the noise 1S more dominant on one phase of the 
line, as shown in Fig. 28, then 60 Hz components of the 
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noise will exist and for a more accurate mode l the 


fundamental period would be 60 Hz. The rate parameter 15s 
Given by 
Ny 
ACC) = E ap coslwo nt) C4 = 3) 
n=9 


where Ny is the order at which the Fourier series is 
truncated. 

The second assumption is arbitrary in terms of the 
distribution function chosen for the impulses. eS 
obvious that random impulse amplitudes are being generated 
by the cumulative effect of the corona noise sources. 
However, it would be difficult to accurately fit a density 
Function to the data due to the time varying nature of the 
process. The independence of the amplitude bursts is also 
conjecture but facilitates the development of a tractable 


model. The Rayleigh density function is 


2 
-a /2B 
pCa) = Ca/Bde a2 @ 
= @ a < @ C4-4) 
and was picked to model the random amplitudes. It 1s 


analytically tractable and can be justified by assuming the 
Quadrature components of the impulse amplitudes are in- 
‘dependent Gaussian events. 

The final assumption 15S made because in most 
interference scenarios, the power line noise process will 


not be synchronized to the process being interfered with. 
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Thus the phase will be effectively randomized over one 
fundamental period. In some of the observations of corona 
noise (see Figs. 26-28) this assumption was clearly violated 
and the instrumentation was intentionally synchronized to 
the power line for clarity of presentation. 

The model could be further generalized by allowing the 
amplitude density function pla) to be a periodic function of 
time, pla,t) independent of the rate parameter. This added 
complexity would make the model more accurate; however, more 
Parameters would have to be estimated. 

The above assumptions classify this process aS a 
compound non-homogeneous filtered Poisson process. A. CGias 
venient tool to study these processes 1S the cumulant 
generating function CRef. 41:p. 117] 


se(t) 
ws) = lntEle 18 (4-5) 


where ln 1S the natural logarithm. The joint cumulant gen- 
erating “Ffunctten.:s 
Ss,e(t,) + Ss2,elta) 
wtS: ,S52) = Int Ele Dio (4-6) 
For the noise process defined by the first two assumptions, 
it 1s shown in Appendix F that 


@ @ 





b ==) Got @, 
ws) = pCa Ata) C e - ljdada 
~@ ~@W 


(4-7) 
and 
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@® @® 
Si ah(t, ~a) + Sea ah(te ~a) 
wCSi ,Se) = pad ACx) C e = ela. 
i be (4-8) 


From the properties of the cumulant generating functions it 


can be shown that 


ECectd)] = w’(@), (4-9) 
VarCe(td] = w’’ CO) (4-10) 
and 
2 
mo wl@,O) 


Cov (t, , te) 
US; USe 


It 1S straightforward to show that 


@ 


Nr 
ECe(t)] = ECal | C an COS (Wo NaIHN(t-adda, 
n=1 


a 


(4-12) 
@ 


a Nr 2 
VarlCe(t)J = Ela J | L~ an COSCWoeNnkIhW (Ct-adda 
n=1 
=m & 
(4-13) 
and 


2 Ney 
SGaovl(t, ,t2ei1 = Ela a | An COS CWe NAINCK, ~-adnNl(ts -adda. 
Nn 


(4-14) 
The autocorrelation Function can be found in terms of the 


above functions as 


RCt, aca = Cov (t, 3 ta ) + ECe C(t, YJIECeCt,)]. €4-15) 
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The above equations, prior to the application of 
assumption 3, define a wide sense cyclostationary process 
CRef. 40] where 

ElLe(t + mTo)] = ELel(t)] (4-16) 
and 
R(t, + mTo ,te, + Milo) = JRCC. ew. (4-17) 
The statistics of the process are invariant to a shift of 
the time origin by integral multiples of the period, To. 
This characterization 1s particularly appropriate when the 
process 1S observed in synchronization with the fundamental 
power line frequency. 

When assumption 3 1S applied, the process start times 

are randomly shifted over one period, the process becomes 


wide sense stationary and the mean and autocorrelation 








become 
To 
a 1 
ECn(t)] = | Sim Cc.) dade (4-18) 
To 
@ 
and 
To 
ee 1 
R(t) = | RiCt+7, tid. (4-19) 
To 
@ 


The power spectral density (PSD) of the process is 


@ 
eS —j wT 
S(€w) = R(tr)e dt. (4-20) 


= @® 
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The following algorithm is an example of how the noise 
process specified above can Pe ditaiiceannlly generated. For the 
algorithm specified below; No is the number of samples in 
one period To, {ni} is the real sequence to be generated, 1 
is the index variable for the sequence and j and k are dummy 
index variables. 

1. 1 = @ 
2. nti) = @ 


3. Generate a Rayleigh random variable, R, with 


parameter AB and set n(i)d = nid) + R 
Ny 
4. Set intensity A = E awzcoslwo ni/No ) 
n=O 


S. Generate a Poisson random variable, P, and set j = 
trunc(P/A + .5) 


6. Set nf€k) = @ for k = 1+ 1 to j 

7.1ies j 

8. Go to Step 3 
The above algorithm will continuously generate a modulated 
Poisson impulse process that may be digitally filtered 
concurrent with the impulse process generation or subsequent 
to the generation of a complete time record. 

Two additional points should be noted. A number of 
impulse amplitudes could stack up on each other if the 
truncation operation returns a zero value. This corresponds 
to the physical case where multiple impulses occur and are 
unresolved within one sampling period. Before a record of 


noise samples generated by the above process is used a 


oA 


truncated uniform random variable between 9 and N,, the sampled 
period, should be subtracted from i. This ensures that the 


start time is uniformly distributed over one period. 


E. PARAMETER ESTIMATION 

Having specified the model, the task of determining the 
parameters for a given interference scenario remains. The 
parameters to be determined are the {a,}, the Fourier Series 
coefficients of the periodically varying rate function and 
the 8 parameter in the probability density function of 
impulse amplitudes. As an example, the parameters for the 
noise process described in the time domain observations 
Section will be determined. 

One method of estimating these parameters is' by 
considering the predicted mean and variance of the process 
as given by Eqns. 4-12 and 4-13 and matching these 
Parameters to the data. In order to do this analysis, the 
corona source has to be observed as a cyclostationary 
process. This requires synchronizing the data collection 
instrumentation to the line frequency so that the same point 
in each cycle can be examined. This is difficult to do for 
field measurements. 

Another option for determining the model parameters is 
the predicted stationary mean, variance and PSD or auto- 
correlation function. The PSD as defined in Eqn. 4-290 is 


stationary so no line synchronization is required and 


22 


Mwerumentation is readily available to make a periodogram 
estimate of the PSD. 
For the measurements taken in this research, the lowpass 
impulse response of the filter was approximately 
2 
ere, 
h(t) = Ca/nie t 2 0 
= @ te GEO. 4-21) 
With the above definition the filter has unity gain at w 
equal to 9. The effective filter response g(t) was approx- 


Imated by h(t). 


In Appendix G it is shown that 


22 
“Wo mn 
1/2 Nr 4a 
Efe(t)] x (Bn/2) Lc hle An COS Cw. nt) A= eo) 
n=9 
ie 
where (nB/2) is the first moment of pla), the probability 


density function of the impulse amplitudes. The stationary 
mean 1S 


ee 
ELeCt)] = (Bn/e2) Ao - €4-23) 


The covariance of e(t) 185 shown to be 
2 
-a(t, -t~e) /e 
Covit; ,te) = (Ca/ndepe 
wo 2 
Ne -2au 
LY an Cos luo NCuU+(t, +t,)/2)e du. 


n=0 
lite = Ce | (4-24) 
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Using Eqns. 4-15 and 4-19 and substituting the above 
results, the stationary autocorrelation uncer enss 
e 


ms -at /2 1/e 2 
R(t) = 2B Ce Gao (a /8n) + (1/4) ao 


ce 
—- Wo Nn 
N+ 4a 2 
+ (1/8) rE e ao Cos Co AT lee 
n=1 
(4-295) 
Taking the Fourier transform, tee resulervaig rsDY i156 
2 
1/e -w /2a 2 
S¢€w) = 28 C2) Sy + (n /@2@daq 6 Cw) 
ewe 
~Wo nN 
N- 4a 2 
+ (€n /8) E e@ Qa S4Cw + nwove 
n=1 
(4-26) 


The terms in the PSD in order are; (1) a white noise term 
multiplied by the frequency response of the Gaussian filter, 
C2) an aiumpulse at w = @ from the mean value of the process 
anes) periodic components at harmonics of 180 Hz due to 
the modulation of the driving fumericon-. 

To estimate the model parameters, the sample mean and 
mean square estimates were used 1n conjunction with Eqns. 4- 
23 and 4-24 to solve for ao. and BB. The periodogram estimate 
of the spectrum was used to determine the magnitude of the 


a; ana a2 co-efficients. The Fourier series was truncated 
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at Ny = 2 since no Aigher order terms were odserved 1n Fig. 
25>. For the data taken in Fig. 24 and 25 the parameter a in 
Eqn. 4-21 was 7.11x1@*%. Table IV lists the parameters of the 


assumed model. 


Coefficient of Gaussian Filter a Toni ee). 
Amplitude Density Parameter 33 aoe los. 
Fourier Co-efficients of Rate es Zens O00 
munCctLlon 
Ay 5, 702 
Aa 2,209 


TABLE IV 


These parameters were inserted into the algorithm described 
previously for generating corona noise, and five cycles of 
the synthetic noise are shown in Fig. 29. Comparing the 
actual noise in Fig. 24 and the generated noise in Fiq. 29, 
the generated noise correctly models the random, impulsive 
character of the actual noise and also 1ts periodicity. The 
synthetic noise, however, does not incor-porate the constant 


Eras above the © level seen in the actual noise. 


F. SUMMARY 

This chapter has presented a detailed analysis of the 
time- and frequency-domain behavior of corona noise. A 
model for synthetic corona no1se was postulated based ona 


filtered non-homogeneous Poisson process. A methodology for 
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Figure 29. Synthetic Corona 


estimating the model parameters was developed. The ad- 
vantage of this model over previous work 1s that 1% allows 
an accurate representation of the time-domain behavior of 
corona. The disadvantage of this model is that it is 
scenario dependent and for each interference situation a new 
set of model parameters must be Re ed. One limitation 
was introduced into the model by the assumption leading to 
Ean. 442. The lowpass equivalent filter is a function of 


the bandpass filter and the process intensity, and 


26 


consequently, the model is not as general or as accurate as 


the gap noise model. 
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ae INTRODUCTION 

Energy detection receivers in the high frequency (CHF) 
radio band operate 1n an impulsive radio noise environment. 
To predict and improve the performance of energy detection 
receivers, accurate models of the signal and noise envir- 
onment must ae used. The noise models developed in earlier 


Chapters will be used to Simulate atmospheric and man-made 


radio noise environments. Two types of energy detection re- 
ceivers will be considered: a fast Fourier transform CFFT) 
processing system anda compressive receiver. TO imMpmeve 


the performance of these receivers the concept Of eeae 
locally optimum receiver 156 introduced to suggest a robust 


post-detection processor. 


Bi» CENERGY DETECTION 
The energy detector recelver measures the energy 1na 
Signal over a specific time interval. For this chapter we 
will consider the signal to be a modulated sinusoid, of un- 
known frequency at the receiver. The signal 1s then 
represented as 
s(t) = Acos(wet + 6(t)) (3-12) 
where A and We. are random variables and #86(t) 15 a Slowly 
varying function so that the signal power 1s effectively 


confined to a narrow bandwidth B. 


he! 


Defining the problem as a hypothesis test 

Flo : a Ge wa Ge wl. noise alone 

ei, : Ce) =e iC tet SCE) Signal plus noise. C352) 
The well known energy detection receiver CRef. 42) uses ten, 


the received energy as the test statistic such that 





- 
[ 2 
Under H.: teyp = | ieee ep oke 
%,) 
+ 
i’ ] 
Under #H, : Cen = Ener S Cl). ee. (SS 2). 
@ 


For a sampled bandpass process the test statistic 15 


N 2 S 

Ho: te o = £ me. ae tiel.guss J 
1=1 
N 2 a 

Bia. s Geo =e EMeavet Sew + Eieueet ose.) CS] 23 
1=1 

and 
hecome = Meu) © Sec 4 Or Acute. ‘Mer 


and likewise for the quadrature component. 

Fig. 30 shows the block diagram of a receiver that cal- 
culates this statistic. THis receiver will be called a 
Square and sum receiver. The square and sum receiver is the 
optimum receiver for detecting random phase signals In 
Gaussian noise at small signal levels. At large signal to 


noise ratios (SNR) this detector 1s very close to optimum 


Shs, 


Hy 











{Xe t } 
Yes 
{Xs t > 
No 
He 
Figure 30. Square and Sum Receiver 
Cher... Sita.. ae Je. When N = BT and the signal 1S sampled at 


the Nyquist rate, the test statistic is the energy 1n the 
Signal 1na period T. 

When it 1s desired to implement the receiver shown in 
Fig. 30 over a wide frequency band relative to the target 
Signal modulation bandwidth B, a number of options may be 
pursued. The channelized receiver approach in whieh iam 
independent receiver operates at each frequency will often 
be too expensive. Two other possible methods for generating 
the test statistic are; C1) a compressive recelver as 2 
spectrum analyzer (CRef. 43] and (2) a fast Fourier transform 


CFFT) signal processor. 


1Q@ 


It has been shown that an FFT signal processor exactly 
implements the test statistic given in Eqn. 5-3 for the 
assumed signal in a Single FFT filter CRef. 44]. Whether 
Successive samples will be independent or not depends on the 
number of transforms per second relative to the siqnal 
bandwidth B. If the sample interval 1S greater than 1/8 


then successive samples will be independent. 


In the compressive recelver, shown in Fig. 31, a fast 
Sweeping local oscillator linearly scans the band of 
interest. Any narrow bandwidth signal in the band will 


appear at the mixer output as a linear, frequency modulated 
ferly, chirp signal. This chirp signal is passed through a 
weighting filter to reduce sidelobes and into a dispersive 
delay line CDODL)I which 1S matched to the inverse of the 
linear FM sweep. The output of the DDL 1s envelope detected 
and, when referred to the Sweep time, it provides) an 
estimate of the spectrum of the input signal. It 165 easily 
shown that for a Slinusoid, the detector  oucput 1s 
proportional to the signal magnitude. Thus, the compressive 
receiver approximates the square root of the ith value of 
fie test statistic in Eqn. 5-3. The samples will be 


independent if the scan revisit time 1S greater than 1/78. 


eee SQUARE AND SUM RECEIVER 
In this Section the energy detection performance of the 


square and sum receiver, shown 1n Fig. 30, will be examined 


1901 
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as the noise statistics depart from a Gaussian distribution. 
The square and sum receiver, which 15 optimum for small 
Signals 1n Gaussian noise, will be used as the reference 
receiver so that the performance of improved receiver 
structures may be evaluated. 

The number of samples, N, chosen for this study was 190. 
This 1S a small Sample size and, for impulsive noise, the 
resulting test statistic will have a distinctly non-Gaussian 
Seescr i bution. In Fig. 32 the probability of false alarm 
(Pez) 15 plotted for unit root-mean-square CRMS) Gaussian 
noise and the Hall models with Vy equal to 2.1, 3, 5 and i2. 
(see Table II). The normalized threshold 1S in reference to 
the unit RMS power of the noise. 

A Monte Carlo simulation technique was used to obtain 
the data for the plot. Since the data are independent 
Bernoulli trials with a large number of trials, the normal 
approximation to the binomial was used to determine the 
confidence interval. For estimates of Pex, greater than 5 x 
POe- with 20,000 trials, the true value of Pe, Wisle! be 
within 20 percent of the estimate with a confidence of $35 
percent CRef. BoD neoe !). The noise samples were generated 
using the inverse method ERef. 495% p:. 951J and the inverse 
cumulative distribution functions C(CDFs) in Table I. 

In Fig. 33 the probability of detection (P,) curves are 
plotted for Pe, equal to 1K1@0-2. A rejection method (CRef. 


45:p. 9521] was used to generate the random deviates since 
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the CDFs of the signal plus noise densities were not found. 
The Pp» is plotted as a function of SNR for the same 
parametric conditions as used in ics Sea. As anticipated 
the Pp decreases dramatically as the noise becomes more 
impulsive. In the case where Vp, equals 12, the Pp decreases 
from greater than 99 percent to 1 percent at a S dB SNR for 
the same Pra. 

A receiver that 1S optimum in Gaussian noise suffers a 
major degradation when the noise statistics become non- 
Gaussian. One possible approach to improving the  per- 
formance 1S to design a rec@éiver that 1s based on the 
statistical parameters of the noise. A parametric approacn 


will be discussed 1n the next section. 


OF LOCALLY OPTIMUM. RECEIVER 

The concept of a locally optimum rec@iver provides a 
general methodology for nonlinear receiver design in the 
presence of non-Gaussian noise. The approach is “Ee 
determine a receiver structure that 1S optimized for the 
difficult small signal cas@ and then examine, and [i956 
necessary, modify the structure for the large siqnal case. 
In this Section the locally optimum receiver for energy 
detection in the presence of Hall type impulsive noise will 
be determined under the assumption of weak signal levels. 
Previous work has considered known signal detection in 


Cauchy noise which 1s related tao one of the Hall models 
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PRetr. 46]. After optimizing the design for weak signal 
levels the performance of the receiver will be examined for 
moderate to large signal levels. The receiver structure 
Will be parametric; that 1s 1% will depend on the noise 
density function parameters. 

To illustrate the method for deriving the locally 
optimum receiver, consider the following Binary hypotheses 
test. We are given a sequence, x1, of N samples. Under «che 


null hypotheses H,, we assume that the sequence consists 


only of independent and identically distributed noise 
observations with a common density pln). Under the 
alternative hypothesis H,, we assume that the sequence 


consists of an additive mixture of a small signal and noise. 
The problem may be stated as 

He: xi = mM noise alone 

rin Kx, = my + AS: signal plus noise Vee pl go N 
with A an unknown small amplitude. 

The optimum receiver, in terms of maximizing the 
probability of detection while keeping the probability of 
false alarm below a certain level 1S given by the Neyman- 
Pearson lemma. This receiver compares the likelihood ratio 
test statistic to a threshold which 1s chosen to achieve the 
desired probability of false alarm. 


The liklihnood ratio test is 


meee x | How > T decide H,, or 
OCA) = (5=4) 
pi ¢x] Hi ) < T decide #H, 
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where p, (x|H. ) ariel Oy ae are the probdability densi, 
Functions for the observation “vector =< sunder ae and en 
respectively. Since the noise samples are independent and 
identically distributed, De (x|Hs > and pryC<|H, 0 “Garnueee 


erry 


written as products of the univariate noise densities: 


N 
Da (x) = Tl one) (S558) 
i=1 
and 
N 
P: Cx) = Il Di Cy = AS, Dig (5-6) 
i=l 
The test can then be written 
N 
HW pi (xX: -- Ase > Hiva@ecide ‘rox, “or 
1=1 
Ax AD = 
N 
i-ey C7 : <t “adeecirde Hye Ge |) 
i=1 
In deriving the locally optimum receiver several 
approaches have been used. Rudnick CRef. 47] derived a form 


of the locally optimum receiver By expanding the likelihood 
ratio in aTaylor series about x and truncating after the 
second order terms. Capon [Ref. 48] showed that the locally 
optimum detector maximizes the slope of probability of 
detection versus signal strength function at the origin and 


employs a statistic based on the following test 


treo (x) = — (5-8) 
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[imeAappendix HH it 1s shown that the racally /) ope rmem 


receiver for the detection of a random sinusoid 15s 


tLe eed = Sur, ? (5-9) 


1 


m2 


1 


where gC:) 1s the zero-memory nonlinearity (2ZMNL) defined 


below 
+ acy) fF’ Cr) 
CC = ee (5-10) 
FCr) rf Cr) 


and r 18S a vector of envelope samples. [Me ecuUnce 1 On. iri) 
for the Hall model of atmospheric noise discussed i1n Chapter 


II is 


(Ce OS gas is xs) 
a = (0-172 
Oe + vy ) 
where K 18S a constant. Prom Sam. o- !O,ene "2ZMNL of the test 


statistic 15S 


e 
C61) CG Soie ato 1) 
Gn). =f SS Ca Yaga 
es ee as eS 
Ge tay) Gan + ¥) 


The receiver that implements this test 15 shown in block 
diagram form in Fig. 34. The 2MNL’s for the unit root-mean- 
square Hall models (see Table II) parameterized by Vp, are 
shown in Figs. 35 and 36. The 2MNL receiver characteristic 
for values of Vp equal 2.1, 3, 6, 19 and 14 are plotted. 

Recalling that the nonlinearity 1S optimized for van- 


ishingly small signals the 2MNL characteristic may be 
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Figure 34. Locally Optimum Receiver for Energy Detection 


intuitively analyzed. Small values of the sampled envelope 
will be given a negative weighting, intermediate values are 
emphasized and fare values of the sampled envelope are 
effectively nulled. As the noise becomes more impulsive, 
which corresponds to a larger value of V5, and smaller @, the 
above features become more pronounced. This 165 seen in Fig. 
35 where the 2MNL has large negative values that extend off 
the plot for very small envelope values. The maximum value 
of the 2MNL occurs at 


Wie 
x Save (6S) / (6 -eie (3 = ie 


and the zero crossing at 


ie 
<=) ¥ (e7 Ce 2) : CS > Ta 


Both of the points are proportional to the mean value of the 


underlying distribution for the same Hall model number. 
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Zero Memory Nonlinearity (2Z2MNL) for Locally Optimum 


Receiver in Hall Noise for Vp 


Figure 35. 
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Zero Memory Nonlinearity (Z2MNL) for Locally Optimum 


Figure 36. 
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=e.1, 


Receiver in Hall Noise for Vp, 


mae sY 1S the Pea versus threshold for the locals, 
optimum receiver shown in Fig. 34. The number N of samples 
1s equal to 10. The Hall noise samples are parameterized 
with the same Vp values used in Fig. 32. The threshold data 


obtained from this plot was used to generate the Pp, curves 


Smownein Fig. 3s. The only complete data set was obtained 
for the Halle case with Vy Squaw eo a. 1. Some data was 
obtained for Vy, equal to 5. The peaks in Pp, for the more 


impulsive noise occur at smaller SNR levels, not shown on 
cme O1lot. The data points are plotted on the graph since 
the dashed curves only approximate the true values. 

It can be seen that the locally optimum receiver does 
improve the P, (Cfor the same Pre,z) in the small signal case. 
mamecomparing the values with Fig. 33 for Vo equal to 2.1 the 
locally optimum receiver has’ a performance better than the 
square and sum receiver from -1@ dB to 7 dB SNR. Above 7 dB 
SNR, the Pp, for the locally optimum receiver drops to @. 
This 1S due to the fact that the receiver is optimized for 
small signals. The Hall2 2MNL in Fig. 36 shows the reason 
mer this. It peaks at approximately 1.3 times the RMS value 
of the noise and approaches @ as the input goes to infinity. 
This behavior of the 2@MNL must be modified to obtain a 
receiver that will also work in the large signal case. In 
the next Section a robust receiver structure 1S proposed 
that does this. A robust receiver 1s differentiated from a 


parametric receiver in that the robust receiver does not 
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Figure 38. P, for Locally Optimum Receiver in Hall Noise as 


Function of Vp for Pez = 1 x 1072 


irS 


depend on the parameters of the noise density and it per- 
forms well, if suboptimally, over a class or even differing 


classes of noise densities. 


E. ADAPTIVE LIMITER 

The 2MNL derived in the previous Section 16 a relatively 
complex funceioen. Furthermore, it does not work for large 
Signal levels. In this section we consider a practical 
implementation of the 2ZMNL consisting of an adaptive limiter 
to approximate the 2MNL function with a modification to 
improve large signal performance. Although the locally 
Optimum receiver was derived for Hall type impulsive noise, 
if the approximation to ic 1S properly designed, 1t should 
show robust performance in the presence of many impulsive 
noise process. 

Examining the Pp, curves of the locally optimum receiver 
for the Hall2 model shows that it provides increased 
detection performance over the square and sum receiver from 
-1@ dB to 7 daB SNR. Since this region 1S where the square 
and sum receiver suffered the worst performance degradation 
in impulsive noise, the Halle 2MNL was chosen to be modified 
to obtain a robust nonlinearity. The solid curve in FPigeges 
is the Hall2 2MNL. To correct its performance at large 
Signal levels, the ZMNL was first altered, as shown by the 
dashed line in Fig. 39, by extending its characteristic 


straight out from 1tS maximum value. Furthermore, an 


Bemrvalent testestatistic is obtained if g’Cx) its any linear 
mmmetion of g(x) CRef. 46] so 


Oo = aac eB a & @. (S215) 
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Figure 39. Modification to ZMNL 


With suitable parameters a and 8, the modified 2@MNL 1s shown 
im Fig. 40. The level, a, remains to be determined. 
Recalling Eqn. 5-13, the argument of the @MNL char- 
acteristic at its maximum value 1S proportional to the mean 
value of the Hall distribution for which 1t was derived (see 
Table I). For the Hall@ distribution the proportionality 
constant is approximately 1.6. Therefore, in order to 
preserve the improvement in Py 1n the same SNR region as the 


Hall2@ 2Z2MNL, the limiter level was set to equal 1.6 times the 


ly 
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Figure 40. Linearized Modification to Z2MNL 


mean value of the assumed noise distribution. 


A block diagram of the receiver is shown in Fig. 41. 
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Figure 41. Receiver Structure for Adaptive Limiter 


Heme cmat the CMNL must be adjusted. This information could 
PeeemiOWN a priori or, as symbolically shown in Fig. 41, 
determined from the eae of the noise. Both cases 
will be examined below. 

The false alarm rate performance of this modified 
receiver in Hall and Gaussian noise, where the statistics 
are known a priori, 1s shown in Fig. 42 for N equal to 10. 
The value for a@ was determined from the noise mean. At a 
Pea @EQGQual to 1 »@ 10°72, the P,p» curves of this recelver are 
Sigowr in Fig. 43. The threshold was determined from Fig. 
42. The Gaussian noise performance of this recelver 156 
poorer than that of the square and sum receiver but only by 
a few percent. However, there 1S a dramatic improvement in 
the Po for the ‘pases hers noise cases. ioe ne sc Time co 
note that for the same noise power, a well designed 
nonlinear receiver provides better detection performance in 
impulsive noise than the optimum receiver for Gaussian noise 
CRef. 49]. 

By plotting the empirically determined threshold from 
Fig. 42 versus the mean of the noise distribution it 1s seen 
that the relationship 1s almost linear. This 1s shown in 
Fig. 44 for two false alarm rates. Using this relationship 
the 2MNL and threshold can be adaptively modified to provide 
nearly constant false alarm rate (CFAR) performance under 
varying noise conditions. This 1S a desirable charac- 


teristic for an energy detection receiver. 


iS 


Pra 


4 


= 


2 


10> ‘ EEE ERE 


TRC 


i nn ee ee eee ee ea 
eS oe 


Ss oe AG ne me 
ale 





Pe ee ee Ee 
a 


See 


mt 


CCCCCERC EN 


| ee 





| 
| 
| 
| 
| 


| 


Tn De ae 
pT aaa ee eee eee 


Figure 4e. 


5 Gq 7 8 
Normalized Threshold 


Pe, for Limiter Receiver in Gauss and Hall Noise 


120 


| 


0. $$39 

j | 
0. $397 | | | 
9. $350 SS ala | / | 


c. S000 / oe). 
0. Sc06 / / J “Mo = 1.5 
Qscca4 ei Vi 


c.3sco4 / Vi 


DETECT TONE PRORAE UE TTY 


C. coO03 
0. cool 
=10..08 ga. 00 ae ale Sy 0 nee a8 PS. 
See ne ee NGiSse KATICG (as) 


cn 
C) 
© 


Eigure 43. Pp, for Limiter Receiver in Gauss and Hall Noise as 
Function of Vp, for Pe, = 1 x 1072 


ai 


= 


-_— —- 





{———- --- — 
a 


-- 


eae eee 


- wee 
- ee 
~ oe 


-_-— 


a 
—_— ee 


-—— -— 


ee 
awe ee 


= 


a ee eee ee 
RE PO es 
a  - — -- 


a fe a ee 


NSS SS 


ee me we Be ee fe 
a = bm 


— eee 
— ee 


te 3 ere | 


= a 


2-97 <X JT pue 
SUOTINGTI4ARSTG TTPH JO ANTeA UPSW SNSUSA PTOUsSeuYy “Hp euNnBry 


1-@O@T X fT 














= ‘Y4q PUP UBATSIOSY UsQtIwiy 


SeNTeA uUedayp 





PToyseuys 











Lee 


A final Question to be resolved in implementing the 
receiver 1s the real time determination of an estimate of 
the mean and mean Square values of the noise. One approach 
1s to sample the Pee in a Signal free adjacent channel and 
another 1s to form the sample mean of the noise using past 
sample values with no signal present (CRef. 5@]. We will 
implement the latter approach and test the adaptive receiver 
on actual: signal and noise samples. In sequencing the 
samples two general approaches may be used. The samples may 
be processed in blocks with no overlap which 1s the case 
corresponding to the hypothesis tests discussed ta this 
Chapter. The other option 1s sequentially processing each 
sample with rank two updates of the test statistic, the 2MNL 
and the threshold. In this case the hypothesis test 
formulation will provide a lower bound on P,. 

The sequential sampling approach was chosen and two runs 
of the adaptive receiver are shown in Figs. 45 and 46. The 
test statistic, shown by the dotted line, was formed from 
the 1@ previous samples using the recelver structure in Fig. 
41. The adaptive threshold, shown by the dashed line, was 
formed from the sample mean and mean square estimates of the 
previous 20 though 5@ samples. The sample mean 15 shown by 
the heavy solid line and the individual envelope data points 
by the light solid line. 

This receiver structure requires HH, to be true at the 


start of the test and will detect the transition to ry 
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Since the adaptive threshold assumes noise only samples, it 
will be corrupted when signal plus noise samples are used 1S. 
FOr Marlies This will occur 240 samples after H, becomes true. 

This effect can be seen in Fig. 45 where an actual 
signal from one filter of a compressive receiver is plotted. 
The noise 1s highly impulsive and would be a very cdi? fae 
detection environment for a receiver optimized for Gaussian 
noise. Samples 5@ to 110 are noise only and samples 119 to 
2e0@ are signal plus noise. The test statistic exceeds the 
threshold at sample 125 and the receiver would declare H, 
true at this pewnc. However, due to the fact that the test 
statistic and the threshold are both corrupted by the signal 
plus noise samples, the test statistic later drops below the 
new threshold. 7 

In ee 46, the receiver performance in a simulated gap 
noise environment 1S analyzed. The signal commenced at sam- 
ple 110 and the receiver declared the signal present at 
sample 11/7. Again the threshold is later corrupted and the 


test statistic drops below the threshold. 


F. SUMMARY 

This Chapter has presented a systematic development of 
the bandpass energy detection problem in impulsive noise. 
Starting with a Hall density function for the noise, the 
locally optimum receiver is’ derived. This receiver 15 


modified to design a practical receiver with robust 
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performance in non-Gaussian noise. The design methodology 
elliearly Hen yeavee the role of the test statistic, threshold 
ana nonlinearity. This allows independent evaluation of 
each of these features. One aspect of energy detection 
receivers not considered is the role of data in filters 
adjacent to the filter being processed. In an impulsive 
noise environment, the noise statistics of adjacent filters 
Will be correlated and their use may potentially offer a 


further performance gain. 


Ley 


VI. CONCLUSION 
Aa RE SUES 

The dissertation has investigated the modeling of 
atmospheric and man-made radio noise and applied the results 
to the bandpass energy detection problem. 

A generalized means of describing an impulsive noise 
process in “termS of its™ complex "envelope "Was derives 
Previous models of atmospheric noise were surveyed and some 
extensions to the Hall model for atmospheric noise were 
developed. 

Using the physical characteristics of gap noise, a nine 
parameter model was constructed that allows arbitrary sour- 
ces of gap noise to be synthesized. The driving impulses 
for the gap noise were postulated to be points of a 
branching renewal process. The application of this type of 
probabilistic model to power line noise 1S unique. 
Additional results include a new proof of the power spectrum 
of an equilibrium renewal point process and a derivation of 


the power spectrum of a truncated branching renewal point 


process. The postulated model was shown to correctly 
predict all the significant features of the observed 
spectrum. 


In a Similar fashion, corona noise was postulated as a 
Filtered impulse process with a periodically modulated rate 


Function. This model was shown to correctly duplicate the 
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observed time domain behavior of corona noise. A method for 
estimating the aaeeineters or the model was discussed and an 
actual example of corona noise simulated. 

The energy detection problem was also considered and 
starting with an assumed density for impulsive noise, a 
design methodology for a parametric and robust receiver was 
developed. This receiver was then tested agalnst actual and 
Simulated data and shown to be superior to the Gaussian 
noise receiver when corrupted with impulsive atmospheric or 
power line noise. 

It has been stated that in engineering design one seeks 
noe, sO much to be optimum but to avoid crippling non- 
optimalities CRef. Sil. It 15 hoped that this research will 
allow system designers to test their systems with simulated 
interference based on the models presented. The robustness 
of the system 1n an actual noise environment can then be 
evaluated, as was done i1n this research for the locally 
optimum energy detector. The noise models developed here 
Should be particularly applicable to systems that operate 
from a fixed site within line of sight of power lines where 
the chance of having interference from gap or corona noise 


sources 1S significant. 


eee FURTARER RESEARCH 
At the conclusion of this dissertation several problems 


are worthy of further study. Research should be done to 


bed 


define the range of values for the power line noise model 
parameters and to determine the effect of power line 
construction practices on the parameters. The suscept- 
1bility of various types of communicatlon equipment to 
Simulated gap and corona noise should also be examined. The 
corona noise model needs to be extended so that 1t can be 
specified and estimated using a bandpass vice lowpass filter 
impulse response. The work on energy detection suggests 
many opportunities for further research. The problem of 
estimating the unknown noise parameters for the adaptive 
updates and the inclusion of adjacent filter samples into an 


algorithm are prime examples. 
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APPEND IXGCA 


INSTRUMENTATION 


The instrumentation configuration employed to provide 
data on the detailed time- and frequency-domain properties 


of high frequency CHF) radio noise 1S Shown in Fig. 47. The 
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Figure 47. Simplified Block Diagram of Measurement System 


radio frequency (RF) input to the instrumentation was either 
a fixed HF long wire antenna or a whip antenna mounted on a 
mobile van. One complete measurement system 1S trans- 
portable and may be moved from site to site. The other 
system 1S installed in a mobile van for noise measurements 
directly at the interference source (Ref. 53]. 

In the lower half of Fig. 47, a Hewlett-Packard 141T 
Spectrum Analyzer 15 used aS a SCaNNing receiver to drive a 


Develco 7200B 3- Axis Display. The spectrum analyzer can be 


jee 


tuned to any desired frequency in the HF = band. Its sean 
rate, scan width, intermediate frequency CIF) Bandwidth, IF ° 
gain, RF attenuation, and other controls can be adjusted to 
best describe the noise under observation. 

An alternate and complementary part of the system 15 
used to examine the narrowband properties of HF noise and 15s 
depicted in the upper half of Fig. 47. An HF receiver 15 
used as an amplifier/translator and tuned to a Frequency 
where the noise is present. The demodulated audio output 1s 
then applied to a Wavetek (Nicolet) UAS@@A spectrum analyzer 
which subdivides the audio-output spectrum into about 5e@ 
frequency segments. When the full audio-output bandwidth of 
S kHz 1S examined, the UAS@@A provides a frequency 
resolution of 1@H2 per segment. The Wavetek analyzer can 
provide individual transforms or average a selected number 
of successive transforms. The transforms Ceither individual 
or averaged) are then presented on the S3-ax1s display or an 
oscilloscope. 

The S-ax1s display provides a continuous moving real- 
time visual representation of the analog output from the 
spectrum analyzers or narrowband receiver audio output. The 
receiver or analyzer output analog data are digitized in the 
display and stored 1n a semiconductor memory. The data in 
the memory are formatted and shown ona cathode ray tube 


CCR) Im aA convenient frequency-amplitude-time (S415) 


ie 


format. The 3-ax1s5 presentation can be frozen at any desired 
time and photographed With a standard oscilloscope camera. 
Fig. 48 illustrates the procedure used to make the 3- 


ax1sS presentation. The analog input 1s divided into Sle 
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Figure 48. Diagram of Data Format 


equally spaced data points Cindicated by the horizontal dots 
or Fig. 48). The signal amplitude at each data point i156 
represented by an 8-bit word. When a scan 1s completed, i1ts 
data are also stored in memory. Line 1 in the view moves to 
line 2, and the new scan appears as line 1. Subsequent 
scans move earlier data, line by line, upwards along the 
time ax1s to create a rising raster type display. When the 
memory 1s full (€64 scans), each new scan 1s entered into the 
bottom line, and the oldest data at the top line 1s 
discarded. The resulting animated view provides a visual 
Picture of no1lse and signals within the block of frequencies 


lf the spectrum analyzer output 1s being observed. becune 
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receiver audio output 1S being observed the display 1s a 
stacked series of consecutive time records. 

The 3-ax1is display has a number of controls to aid the 
observer in analyzing the Structure of data presented. 
Among these controls are: (1) elevation and azimuth geometry 
controls to vary the viewing aspect, C2) an amplitude 
compression control, (3) a threshold control to vane 
background levels, (4) time-axis expansion controls, and (5) 
a stop-action switch to freeze the data in memory for 
detailed observation, For observation from various aspects, 
or for ~phoGegrearpihine . These controls are used to optimize 
the presentation of desired signal detail. 

Accurate frequency, time and amplitude calibrations are 
maintained so that the resulting 3-axis photos can be 
manually scaled for precise signal detail. The digitized 


data at the display input are available at a digital 


interface connector for external digital recording #62 
processing. Received data can be processed By computer to 
obtain conventional noise statistics. This feature permits 


the comparison of selected 3-ax1s vilews of noise with 


standard statistical descriptors of noise. 
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APPENDIX B 


So TERED IMPUBSE PROCESS 


The complex envelope of a noise process can be described 


by the following equation 


Not) aefeln 
net) = g(t) + XL age Ce © aD (Bat) 
1=1 
where g(t) 1S aA Quadrature Gaussian process, Bee, 1s a 


linear filter, a: 15 an independent, identically distributed 
random amplitude and 6, 1S aA random phase uniformly 
distributed over @ to @n. The joint characteristic function 
of the inphase and quadrature components of the above 
process 15S 


jw: Me CO) + Wats C(O) ) 
DCw: ,Wwe) = Ele le (S—2) 


If the Gaussian noise 1S assumed independent of the 
impulsive noise the characteristic function of the Gaussian 
term alone can be written 


3 OW: Oe Ce.) + We Os Ct) ) 
Da CW 9 We ) = ECe ale R= so. 


The inphase and quadrature components of bandpass white 


Gaussian noise are independent so 
a 2 
-N, BCuw, + We ,/2 
Da (Wi ,We) = @ ‘ (B-4) 


The joint characteristic function of the ith pulse of the 


inphase and quadrature impulsive noise terms i15 


1 eS) 


jar (We Ge (Eb -Ci J ECGSO7) sahaeCe Cpe Sc ee 
OQ, Cui ,we) = Ele 


Wa C Re Cet, Is1nNO, + Rs CG—c, Ycosé, )) 


lg: (Bl =p 
Since the noise process is narrowband, its” joint densre, 
funeeron' 1s circularly symmetric and the Following 
transformation can be applied 
eo fan Oe Shas sal 
Wr = Cu, + we ) a = tan Cw:i /Wa ) 
WY: = WwW, COSa Wa = W,S1INa. CBS? 


This type of transformation of a 2-D Fourier transform to a 
one dimensional transform of the envelope is known as a 
Hankel Transform and the inversion formula 15 
ey 
ptr) = | ru, OCO, J Jo (rw. ae (Baa 
@ 
where Jo is the ordinary Bessel function of the first kind, 
order @. To simplify the derivation it 1s assumed that 
ankte ) 1s zero (see comment before Eqn. lies 1) Now 
transforming Eqn. B-S and taking the expected value of 8, 
and a, 
o an 


f i jWn ar Ne CE=; 2Cos CO a) 
Q@, (w,.) = E¢| [Pca pCa re dé@,da, ] 
Q@ @ (Be) 


Since 0, is uniformly distriGuted over oven, 


eer 


Qiao) —se twa Cay G6 (0, a, Me (7) Ida; J]. 


Q 


CB-—S) 


re MCT) in Eqn. B-1 165 a Poisson process with rate A and 


observation interval T, the ™ ’S are independent and 


over (@, T) allowing the characteristic function 
impulsive term to be written as 
4) K 


O0,¢€wr.) = £F OBO Cw.) PrtNnN = kiI. 
k=@ 


Now using the probability law for k events of a 


process 


foe) kK -AT kK 
2,(€wre) = EX BO Cw.) e CAT) /k! 
k=@ 


which sums to 
@ Wy 


A) Pea» | 55 Coe abe 9) - Ljidtda 


@ @ 
Oy, Cw, ) = e ° 


Using the convolution property Or = zero order 


uniform 


of the 


(B=) 


Poisson 


CE=t) 


CB=1 2) 


Hankel 


transforms CRef. S41] and Eqn. B-7 the density function of 


the envelope for Gaussian plus impulsive noise 15 
@ 
ptr) = rw, Jo (rwr )GDe (wr ID: Cw, Ddw,. 


™~& 


7, 


CB=13) 


APPENDIX C 
SIGNAL PLUS NOISE Seehs Ty 
For the truncated Hall model with 8 = 3, the joint 
density function of the in-phase and quadrature components 


1S given by: 


2 a 
D - 2y a 2 2 
p(x,y) = Oo > Se 
Fes 2 2 Be 
en(D =1) (x ya 
es 2 = 
= @ x Hye > wie (Cl 
where 
Flere 1G) 6 Fa 


Dea CPS ey yo) : 
The sum of a random phase sinusoidal signal of amplitude A 
and noise | can be expressed in terms of its in-phase and 


Quadrature components w and z where: 


Ww Acos(w) + x 


z = AsSin(w) + y. (C=2) 


Eqn. C-2 is substituted into Eqn. C-1 and transformedies 
polar coordinates with 

w = rcos(d) 

Z = rsin(O) (C =p 
where r represents the envelope and 2 the phase of the 
Signal plus noise. The joint probability density function 


p(r,O) conditioned on w is 
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B 2 
D rey 


p(r,@| yw 
e 2 2 etic 
en Gem) (ice —-eeArcos Cie - DO) + Att y ) 


Q a a 
Oo < riiemaecArcos(w - 0) + An oc Te 


= @ elsewhere. (C=) 
This expression is valid for all values of r between @ and 
iene - Az Since our principal interest is the small signal 
case, A << Te, this will not affect subsequent analysis. By 
integrating Eqn. C-4 with respect to @, the conditional mar- 
gQinal density function of the envelope, p(r{w), 15 obtained; 


n 


“\ 


PCr) P= se) por .? lw) doe CG=5) 

@ 
The integral is a function only of costw - 8), w appears 
nowhere else in the integral. Therefore by setting g0 = wy - 


®@ and using Eqn. 2.5594-3 in CRef. 54] gives 


2 a 2 2 <2 
Cir +e pete y Daev: iD 
p(riy) = 


4 22 eae (ad Ws 4 4 2 
COmaemhe ten + CA Y +eeyv rm +t At y CD = 1) 
CE-—o> 
Since w 1S assumed independent from sample to sample and 
uniformly distributed over en, 
en 
p(r) = | p(r{[wptw) dy (C-7) 


@ 
and 


ESSE, 


zz 2 2 Ee ie 
(+ A Geeave) 2 yaa 
2 a) 
4 ema ae ee £ 4 2 
Cr - CA fm + GA VY +° Gyr ete ye Cee 
(C=e) 
By a Similar argument, the density function for a random 


phase sinusoidal signal i1n narrowband noise with a Halls 


distribution 15S 


4 2 2 
ey r(G + 2am 
1) SS E's) ) 
2 2 ere 
Ca le) 
where 
a 2 a 
a= r +A + Y 
and 
b= =-=2Ar. 
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ARFrEND 0 7 


SEFeciaRum OF A RENEWAL PROCESS 


The complex envelope representation of filtered impulse 
noise ona finite interval T 1S given by 
NCT) 0% 
n(t) = LF ae Cesc lus Daa” ASB. 
1=1 
GO = 1.) 


NCT) 1S a random unit counting process denoting the number 


of impulses in the interval, the {t,} are the random arrival 


times and 98, is uniformly distributed over @ to en. The 
pulse amplitudes Ca, } are identically distributed 
statistically independent random variables with second 


moment Az and fourth moment A,. The filtering effect on the 
impulses are accounted for by the time invariant complex 
envelope impulse response H(t). 
The envelope squared of n(t) 15 
2 NiCL) ce 
EGcoee— MCG Mm Cll. = se a, NCe-t, Ihc -t, 
i=1 
CD32) 

where we have assumed terms of the form HACt-t, Dh" (t-t.)9 for 
different i and kK are negligible. We estimate the power 
spectral density CESBo= “oF the envelope By averaging 


magnitude squared, lenge i.) -OUrieretiransrorms of E* (t). 


The mean value of the estimate 15 


Sti) — ELS. CoS. co) J 771 (Pe S)) 
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where 


+ 
; 2 = joc N oT) 2 = tor 
1=]1 
Q 
and -@ 
=e 
He (wd = nce lh" Cee ac 
@ 


1s the Fourier transform of the magnitude squared impulse 
response. We have assumed the impulse response 1S much 
shorter than the observation interval T. 
The mean of the estimate becomes 
_ NOCTIUNCT) “ae.2 SiGe ree (i 
S Coe) = ise es E-Ga; au Ue Vine eooy 712 
k=1 i=1 | 
(D-4) 
This 1S a compact and general expression for the PSD of a 
truncated filtered point process. It 15 1n terms of a 
summation a al the characteristic Functions Ore the 
interarrival time distributions between all combinations of 
pairs of points. (The term incorporating the frequency 
response of the filter will be set to one for the rest of 
the analysis.) Conditioning on NCT) 
2. Creating. — Cy) 


N 
oe QA: a. e IPriN ="Tivrte 


on se) N 
pate Cited BD tee — es ae Em 
= = | 


Te 
(D=27 


Plotting the difference k-j indicates a different summation 


over m and n shown below. 
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Let m = 1 Stee ang ny — lett such Chat ct. = t, - tt. is the 
mth interarrival time. Then Eqn. D-S can be written as 
_ rs) Ned, =n ee wee 
Stw,l) = <£ Ef ZL x“ Aa e + NA, 
l=2 n=1 m=1 
eee en 4 te 
os x“ As e PRrmEN = iI7 i. 
n=1 m=1 (Db-6) 


On further rearrangement 
= = oon N-1 n jwt. a ies 
SCw,fI= MA,/F + E£ Asp EL £ roe ae = dae Nee bd / ile 
l=2 n=1i m=1 
Coe) 
and conditioning the interarrival time distribution over the 
number of pulses N in the interval T 
= = Co LS n 
S€w,T) = NA,/T + Asa CL = © CO. Cjw|/N=1) + 
l=2 n=1 m=1 
Oo." (jw|N=1) IPrtCN = 1I/T 


(D-8) 
where 


SUS 
@.(j;w}/N=1) = ECe } 


is the characteristic function of the mt" interpulse spacing 
Qiven exactly 1 pulses inthe interval. The following 


assumption 1S now made; that the mth interpulse spacing is 
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the sum of m statistically independent interpulse  sep- 
arations, and 1S independent of 1, the number of pulses in 
the interval. TAisS assumption will give an asymptotically 
correct result for the PSD as the observation interval T, 
becomes much greater than wp, the mean value of the 
interarrival times. 

With the above assumption 


mM 
De Co | N= Tee Cee sD (D-9) 


where @( jw) 156 the characteristic function of the iunter- 


pulse separations. Eqn. D-8S then sums to 
n+} 
= = 2eo 1-1 @O(j0)-2 Gye 
Stw,T) = NAg/T + Aa £ — 
l=2 n=1 1-O¢€ jw) 
na 


OFC 3w)-O" ( joo 
ss EEN = ee 


1=O8°C wD 
(D-10) 
and summing again 
| 
= = (fake es) 1d © 5) QC jw) OC jw) 
S€w, 7) = NAg /T + Aa Bo 0m se eee 
l=2 a 2 
C1-O¢ jw) ) Ci=2 Ca C1-OC jw)) 
oe 
1@" (jw) @" ( jw) @" €3w) 
7 Sree +Prt(N = Pave 
a a 
Cho" 2G io) C1=0”" (ao) Gi OraGaa) ) 
(b= 1) 
Rearranging 
= = 2 = ee) OC jw) @" € jw) 
S€w,T) = NAs /T + Aa STAN = Prt N=1i) 3. —— See) coca 
l=2 


(1-O(jw)) (1-0(€5w)) 
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a OC jw) DB" (jw) 


ee See Ser 2) le Nl a ——— SCO 
@ Va 
C=O Cw) GE=Os C ja)? 
Let Ded 
zs ee) QC jw) O" (jw) 
toa /) = Prt{N = 1] rr Cis 
l=e 2 ze 

Cr -O Gu) C—O“ 7G aD GCD= 12) 

Since Prf{N = 1] is less than 1 and {@(w)}] < 1 CRef. 


ae Vio form amy distribution, then the absolute 
convergence of the third term 1S guaranteed. For an 
equilibrium renewal process 

N = ECN(t)] = T/p Coen) 
where p 1S the mean time between renewals CRef. 3e:p. 46]. 
Now letting 7T, the observation interval, go to infinity the 


asymptotic PSD 1s 





a = Aa As DC jw) 
S€w) = lim Sw,T) = ea tL ee 
T=>0 y Z 1-O¢ jw) 
Aa Cis ) 


This 18 the expression for the PSD of an independent 
increment point process in terms of the characteristic 
Function of the interarrival distribution. A similar result 
result was developed by Cox and Miller by considering the 


renewal intensity function (CRef. SS]. 
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APPENDIX €& 


SPECTRUM OF A BRANCHING RENEWAL POINT PReges> 


The complex envelope representation of filtered impulse 
no1se arising from a branching renewal process on a finite 
interval T 15s given by 

MCT) Ne eae 
net) = CF PA ok Oe Ge Steen = a (E> 
m=1 n=1 
where T, 1s the beginning of the mth main interval, Nw 15 
the number of impulses in the mth interval, arriving at 
times {ta,}+} after 1ts onset and MCT) 18s the unit counting 
process defining the number of main events in the 
observation interval TT. The pulse amplitudes {a.,} #£xand 
-phases {(8.,)} are statistically independent random variables 
and the {8.,}) are taken to be uniform on {@,2n}. The net 
filtering effects on the impulses are accounted for by the 
time invariant complex envelope impulse response, NAC(t). 
The envelope squared of the process 15 
2 MCT) Ne 2 (a: 
E (@)-= nce yn ce), == 2 ye ae a eo ee i) | (E-2) 
m=1 n=1 
and using the same assumptions as Appendix D the mean value 


of the Bartlett estimate of the PSD is 


= 2 
S(w,T) = EC{Es(w,T)| I/T (E-3) 


where €2 (w,T) is the length T Fourier Transform of E? (t). 


specifically, 
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None = a Wicars = OES 


with Ha(w) = | OR NY the Fourier transform of the 
@ 
magnitude squared impulse response. To find the mean of the 
estimate we must evaluate 
_ eermiCT = Neo MCh NN, 
Stw,T) = EC |He Cw) re = NE, az 


m=1 n=1 1]=1 k=i 


= ion = {0 Goan =e natasha. D 


Aan Aix e 2 I7 i 
CESS) 
Im erder to simplify this expression, we make use of 
assumptions i, 2and 4 from section 3D3 (1) Ela®?.,a?1,,.] = 


a*, which will be normalized to one for simplicity, (2) N, 
and N, equal either Ni or Ne alternating between the two 
values and GS age (Goa) and f(t,nA), the interarrival time 
distributions, alternate between the negative and positive 
phase distribution parameters. Furthermore, let M be an 
even number of malin process points occurring at intervals, 
To/2, Of a fundamental frequency (Cin our case f. = wo/en = 
1/T. = 6@ Hz, the power-line frequency) and let Me = M/e. 
This assumption will greatly simplify the calculations and 
for M greater than 190, not significantly reduce the accuracy 


of the result. 
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Now looking only at the factors involved in the 


summation in Eqn. E-5S, six separate terms can be identified: 


Term 1 - The inter-pulse group summations between pulses i1n 


pulse groups with Nl pulses 


Ni Ni jC Cae ene 
= Mer = “Efe J 
n=1 k=] 
(E-6a) 
Term 2 - The inter-pulse group summations between pulses in 
pulse groups with Ne pulses 
N2 Ne jw Cont ae 
= M2 Ct ce ce i) 
n=1 k=1 
(E=6ise 
Term 3 - The intra-pulse group summations between pulses in 
groups with Ni pulses 
M2 M2) =jot ma lem. Ni Ni pee =e. 
= -L x oe 1B r Ele JECe J 
m=1 l1=1 n=1 k=1 
mn< | 
(E-6en 
Term 4 - The intra-pulse group summations between pulses in 
groups with Ne pulses 
M2. M2 =—joCGCm=1oi. N2 Ne Hes = Wc. 
= - ee: DE mele JECe J 
m=1 l]=1 n=1 k=1 
m<>1 
(E-Gae 
Term 5S - The intra-pulse group summations between pulses in 


groups with Ni pulses per group to pulses i1n groups with Ne 


pulses 
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Mota = ip ame Me ~-jwCm-1)7T. Ni Ne Fie. aay CO Gy 


=e De x he DE 2 Ete JECe J 
m=1 l1=1 n=1 k=1 

(E-6e) - 

Term 6 - The intra-pulse group Summations between pulses in 


groups with Ne pulses per group to pulses in groups with Nil 


pulses 
macs /e-1p) Me me = Mele) are N2 Nil Lvicn -jwt. 
=e L “oe x zZ ele JECe ] 
m=1 l=1 n=1 k=1 
CEaGt> 
ft 1S important to note the expression for the 


characteristic function of the interarrival times differs 
between the inter-pulse and intra-pulse cases. In the 
inter-pulse case, discussed in Appendix D, the summation 15 
only over a function of the difference in the pulse 
positions where in the intra-pulse case the summation “ 
over the absolute position of each pulse i1n 1ts group. 

Terms 1 ande2 can be evaluated using the results from 


Appendix D and assumption 3 from section 3D: 


Term 1 = Me S(u,N1) CE-Ja) 


Term 2 M2 2¢€w,Ne) CE— 7b) 
and the remaining terms can be determined by straightforward 


evaluation of the finite sums 





GosomMe enw/w,2>- 1 
Per 38) iO OA .——— Pe de APC ye fal, 
cos(2enw/w.) - 1 C2=7e) 
cos(Me2e enw/we) - | 
Term 4 = ¢ - - M2 } Q(w,N2) 
cos(e2nw/w.e) - i CE=7d) 
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Terms S and 6 are complex conjugates and can be combined to 


form one term 


cos(Me enw/we) - 1 
Term Sa = 2cosl(wtl, 42-10) SSS ee Nee 
cosCenw/w.) - 1 
CE-7e) 
where @, and @2 are the characteristic functions associated 


with the interarrival times for each phase and N1 and Ne are 
the number of renewals associated with each phase. The 


functions &, 2, and w are defined below; 


2 N 
N - N{S@¢Cjw) | OC jw Grae Cue 
S00 N= 
2 2 
eee | CA ioD) a=) 
eral Fe 
[BWC jw) - BOC jw) | 
I 
2 
[1 - @C jw) | GB) 
and 
Ni+1 N2+1 
C(O, Cjwd - BO, C jw) ) (Oe - jw) = Os C= oF ) 
wiS, ,N1,o, ,Ne2) = esses ______nniiindi;nnnsinin;,;,:+:iaae 
C1 = @y (C507) Cle oo oe 
CE-1@) 
Thus 
i > 
S(€w,1) = tECHs Cw) 1/7? 2 fermcn (E= i 


n=1 
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APPENDIX F 


NGM HOPICGENEOUS POLSSON PROCESS 


The following derivation is for the mean and variance of 
a compound periodic non-homogeneous filtered Poisson point 
process. hls aroor 1S a generalization of a proof 
presented in Papoulis CRef. 4l:pp. 382-383]. The expression 
for such a process on an interval @ to T 1S given by 
N¢T) 
ae) = ey ieee (F-1) 
i=1 
where NCCT) 1S a unit counting process with periodic rate 
parameter Att), a:°1is the weighting of the ith point and is 
independent and identically distributed CIID) from point to 
point and ACt) is the time invariant impulse response of the 
wLlter. | 
The time axis 1S divided into consecutive intervals I, 
Peerength Aa where am = @, - &,-;1.- The number of jumps in 
the counting process in the jt" interval [, 15S given by am,. 


Pear 1s SUFfFICIently small, then the contribution to the 


total noise process from the jt" interval 15 


am, 
Anyeese = a, PGE-a, ),; GE ) 
1=1 
where am, 1S a Poisson random variable with 1ts rate 


parameter approximately equal to 


ICO, Joe CE = 30 


Lol 


The moment function of the jt" interval is then 


am, 
= fore Baer= WW eh e, 
1=1 
40, (Ss) = Ele ate 


Conditioning on an, , 


AM, 
s kL ai gGaeeaieds ) 
1=1 
4@,(s) = Ele fam, = k] Prfam, = ki] 
and 
kK Sa, ACt-a, ) 
40,(s) = EC MN e ] Priam, = 


L=—1 


Since the a, ’S are independent from pulse to pulse 


kK sah(t-a, ) 
AS, (Ss) = Si Ere 1 Prlam, = ki, 
1=1 
where the expectation 1S with respect to the 
amplitude a. Using the Poisson probability law and 


nizing the series expansion for the exponential gives 


sah(t-a, ) 
ACa, JaxtELe i a> 
49,(S) * e : 
Since 
@ 
NCE) = -L .“AniezGeoe 
1=-o 


a sum of independent random variables, then 


wm 
O®,(s) = IT a®, €s) 


1=-o@ 


(F-4) 


@ te), 


(Fee 


( Fae? > 


random 


recog- 


(FP =s) 


CP <3) 


CF =109 


using the convolution property of the moment generating 


function. The cumulant generating function is defined as 


ie 


eS) eGo. (S)), Cae Vee, 


where ln is the natural logarithm and 





wm 
Wa (S$) = EL InlaO®, (Ss)]. Ceo 12) 
1=-o 
AS aa => @, then 
f+ a] @ 
i sah (t-a) 
Wn CS) = pCa) ACa) Ce =sliJidada. Gr a1'3)) 
~ wD nae ¢ 4) 


Using a similar argument the joint cumulant generating 
Function for two random variables, n(t:) and n(ta) 15 
@w @ 
; S,aNh(t, -a)d) + S2gah(ta -a) 


Wa (S1 ,Se) = } pcta)d | ACa) C e - lidada. 


“@w ~@m 


. (F-14) 
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APPENDIX G 


AUTOCORRELATION OF CORONA NOISE 


The mean value of the cyclostationary process defined by 
assumptions 1 and 2 in section 40 1S given by Eqn. 4-1e. 
Substituting Eqn. 4-20 gives 

ee) 2 


+. | 1/2 -as 


ECe(t)] = ELlald Aan COS CWe NC t-S)) Ca/n) e ds. (G-1) 


ee ad 


n=@ 


@ 


Expanding the cosine and integrating 


ene 
-We Nn 
Ny 4a 
ECettd)J] = (Bn/4) © C e Aan cosl(wWo nt) + 
n=O . 
ae 
Wen 4) 1 -Wo Nn K= 1 
———————— = —S = C0 Aan SINCWe Nt) J 
1/2 k=1 (C2k-1)!! 2a 
2ctcan) 
(Gaee 
72 
where (Bn/2) = ElajJ and ¢(2k-1)!! = 13 5...2k-1 (CRef. 34: 
PGi Se 897 Jes The sine term 15S 1n Quadrature with the 


dominant cosine term of Eqn. G-e and represents a phase 
shift in the mean of the process relative to the rate 
FUnGE lon. This 1S due to the filtering effect where the 
impulse response of the filter persists, thus causing a lag 
in the mean value of the process. For the values of a and 


Wo considered in this work, the sine term in Eqn. G-2 will 
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be negligible relative to the cosine term. Piivone auiey tints 
means that for filter bandwidths much greater than the 
Fundamental frequency of the corona noise, the phase shift 
will be insignificant. 

Similarly the variance of the envelope when considered 


as a cyclostationary process 15S 


Ny 1 ; ie Sa 
Varfle(t)] = 28 DL eS oe ED 2 an eos (Wo mt J. 
n=O = 
CG-3) 


Using Eqn. 4-14 the covariance of the assumed process 15 


o 2 & 
Ny —CCC = S) =“ CEs —S) 
@oyv¢t;, ,te) = 2 | L an coslWo ns) (a/node 2 ds 
Jv n=O j 
max Ct, , ta) CG-4) 


where after rearranging and completing the square 


& 
=~ (CC, -t2) /2 
Cov(t,; ,te) = (2RCa/nde ) x 
~ e 
Ny = Bo (ey (eyes Dye. 
LF 8 an COosS(We.nsdJe ds. 
n=@ 
maxCt, , ta) CES >, 
Letting 
ie Gee CE te) 2). (ESS). 


the covariance i156 


LoS 


(a 
I (Chea ae 


COV Ct, .t22 (Ca/n)2Re > x 
ee Find 
Ny = Ul 
YT an COS Co mCU+-Ce, cy ae du 
vy n=@ 
(t, -ta | CG=7) 


which cannot be evaluated in closed form. 

The stationary autocorrelation function 1s Foltndae 
using Eqns. G-2 and G-?7 and the definitions given by Eqns. 
A= and 4-19; The term due to the covariance in the auto- 


correlation function is 





Fon lee @ 2 
-at /2 C Ny -2au 
CR Ca/ndoe Fate | L£ Ca, COS (We NCL2C+tT) dD Cos (We nude Gia 
J n=O 
Q {7 | 
co 2 
i -2au 
+ a, SlMCWso MC cette? Sin(€Wo nude Gliamiic tr 
as (G ser 


The sinusoidal terms will integrate to zero with the 
exception of the ap cosine term and the expression simp- 


lafies- te 
Fd ee a 
wel 8 in Cae ams = faeuw 
CR Ca/nde Ao e Cher. 
| | 
(G-9) 


This term will be dominated by the first exponential and 15 


approximated by 
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2 
-at /2 1/2 
cCRBe Ao Ca/S8n) . (G-10) 


iie stationary autocorrelation function term due to 
cyclostationary mean value function is 


To 
1 
= | SiGe Gey JECe cert). dt, CGiz tl) 
To 

@ 


which using Eqn. G-2 1s equal to 


22 
-Wo Nn 


T 4a rad 2 


2B (1/78) e Aan COSWos nT + 2RC1/4) a, . 


m2 


CG" 22 


An approximate expression for the stationary autocorrelation 


> 


Function is 


2 
= -at /e 1/2 2 
R(t) = 2BlLe Ao (a/8n) + (1/4) ao 


Nr 4a a 
+  Cl7ee L oe Aan cCOoS(Went)]). 
n=1 
(G-13) 
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APPENDIX H 


LOCALLY OPTIMUM ENERGY REGEI VER 


The observations consist of a sequence of N complex 


samples. The hypotheses are: 


Flos: Ky = 
versus 
Ne ; 
KH, 3: x4 = mri aA A Si 1 = oS SAN C(H-1) 
where m are an independent, identically distributed (IID) 


complex noise samples with a joint inphase and quadrature 


density pl(mne,Mms). The sequence Ss, 15 a complex, zero mean 
Signal sequence with a known variance, a:, and Aisa real, 
positive number. The 1/2 power of A 1s chosen to facilitate 


the derivation. 
The generalized likelihood ratio (CRef. 6:p. 585,21] to 


test the hypotheses described by Eqn. H-1 1s 


fe @ wo 
1/e Pe 
N D(xe: 7A Sei yXsi 7A Sai) 
Atx,A) = P(Sce:i ,S81)40Sei1dSe,. 
1=1 Dxe1 »Xs1i ) 
-~ -@ (H=2) 


Applying Eqn. 5-8, and differentiating with respect to A, 


s t/e ive 
-Se1 ONC ei A Sei 9X01 7A Ssi)/e 


1 ive 
OXe1 DO&e 4 5 <s 4 IA 
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i7 2 i7e2 
—-Sei1 TP Ox%ei-A Sei sXs17A Sai dle 
J plSe1i ,Ss81 2dSe.1 dSsi 
i172 A=>@. 
OXei1POx%e1,Xs1iIA 
Cas) 


Letting A go to @ and simultaneously taking the expectation 
Of Sc. and Ss. requires L’hospital’s rule to evaluate the 
expression [CRef. S56]. After taking derivatives the test 


statistic is 


@ 4) 
\ exe 7 2 L/G 
N Sei 0 Pl®Ke1 -A Sei5Xs1 7A Ssi)/2 
tio (x) =E£ at 
1=1 a 
OXe1 PCXer > Xai) 
-o@ -@ 
2 1/2 lire 
Sei1Ssi10 Plxe:1 “A Sei 5X31 7A Sei) 
SR sf 
OXe1 OXe1 POXe1 5 Xai) 
oe We ive 
Sei: 0 PlX%e1 -A Sei sXs17A Ssi)d/2 
FT née, S01 )d5c1 Ss 
2 A=>@. 


OXsi PC&e:i 5 Xsi ) 
CH-4) 


In general, the inphase and quadrature noise components will 
be circularly symmetric even though not necessarily 
statistically independent (Ref. 41:p. 133] such that 
(i 2.1/2 
f(r) = pl €xe + Xe 0) ) CaS): 
with 


va Venda 1 Aad 
r= (xe + Xs ) < CH-6) 


UR 


Now A may be set to @ and 


@ @ 
2 Ved 2 2 
N Sey Ery f° Cry xen = ee 6 er ee 
teat.) Wk ae C rs 
1=1 3 
ae Ce 
“Oo ~@& 
Se, Sei Cri’ ’ Cr, ee 1 xe OF FO Or ec so 
+ 
3 
ei tra) 
Fa (oa a 2 
Sea Crypt? “Cr, Der = hae Ce on tea ee 
Jp(Se;1 ,S81IdSe, dSo;,. 
3 
re fF Cire CHa 7) 
If Se. and Ss; are assumed to have equal variances and a 


Zero covariance the test statistic 165 now written as 





N 2 
C6 (i) k= D4 (Oe Oe CH=S2 
1=1 
where 
ft ~Gigd) Cie) 
g(r) = - —————— . CHS) 
ne Me, int ia 


If the unknown signal components are assumed to have equal 
variances from sample to sample then an equivalent test 


StacrecLic 1¢ 


© ues Ge: Be © AG aye A CH= Tap 
= 


16@ 
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